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I.  INTRODUCTION 
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The  oceans  of  the  world  are  an  aggressive  environment, 
both  towards  man  and  the  equipment  he  operates  within  them. 
Prediction  of  the  deleterious  effect  of  sea  water  on  metals 
is  but  one  facet  of  man's  attempt  to  adapt  to  that  environment. 
The  practice  of  over-designing  must  end;  we  must  design  to 
legitimate  specifications,  not  over-design  to  mask  areas  of 
ignorance.  Corrosion  is  one  area  where  much  information  is 
yet  needed.  Determining  corrosion  rates  is  a complex 
problem  indeed.  They  are  the  product  of  a myriad  of  para- 
meters, some  controllable  by  man  and  some  not;  the  velocity 
of  the  sea  water  past  the  metal  surface  is  one  parameter 
that  may  be  either. 

The  effect  of  increased  velocity  on  the  corrosion  rate  of 
metals  has  been  studied  for  over  25  years  [l].  In  all  that 
time  the  most  serious  shortcoming  of  such  studies,  whether 
with  rotating  specimens , revolving  electrodes  or  flowing 
electrolyte  systems  has  been  their  inability  to  be  correlated, 
and  the  inability  to  apply  data  from  experiments  to  service 
situations.  Table  1 and  Figure  1 show  the  results  of  work 
by  several  investigators  to  determine  the  relationship 
between  velocJ  cy  and  corrosion  rate . While  the  trend  is 
obvious , so  too  is  the  disparity  between  reported  values  . 

Some  of  the  scatter  may  be  due  to  the  use  of  small 
samples  where  non-uniform  corrosion  (crevice  corrosion, 
leading  edge  attack  or  pitting)  could  yield  skewed  results. 
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Corrosion  rates  reported  for  90/10  Cu-Ni  in  flowing  sea  water. 
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Figure  1.  Variation  of  corrosion  rate  with  sea  water  velocity  as 
found  by  several  investigators. 


But  the  problem  of  correlating  corrosion  rates  in  different 
flow  situations  involves  other  factors,  particularly  those 
associated  with  characterization  of  the  flow  itself,  and  in 
some  cases  these  (flow)  effects  may  dominate  the  corrosion 
behavior,  rather  than  other  factors  associated  with  the 
environment . 

The  character  of  flow  near  the  metal  surface  is  of  prime 
importance  because  it  represents  the  major  difference  between 
various  test  setups  or  service  situations;  for  example,  the 
difference  between  flow  past  a small  water  wheel  sample 
and  that  past  a high  speed  hydrofoil  craft  [8].  Historically, 
the  corrosion  rate  (R)  has  been  plotted  simply  against 
average  fluid  velocity  (u) . Porte,  for  example,  related 
the  two  by  the  equation: 

R = + k2Cu’^ 

where  C is  the  oxygen  concentration,  k^  and  k^  experimen- 
tally determined  constants,  and  the  exponent  n,  in  the  range 
1/3  to  1,  depending  on  the  type  of  flow  and  surface  geometry. 
The  first  term  is  the  natural  convection  contribution  which 
is  insignificant  at  high  velocities  [9], 

To  date,  attempts  at  relating  flow  characterization  to 
corrosion  with  non-dimensional  groupings  of  parameters  have 
stopped  with  the  Reynolds  number 

Re  = 

V 
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where  u is  the  average  velocity,  L is  some  characteristic 
dimension  (distance  in  the  line  of  flow  for  a flat  plate 
and  inside  diameter  for  tube  flow)  and  i'  is  the  kinematic 
viscosity.  Conjecture  on  the  use  of  other  groups,  both 
dimensional  and  non-dimensional,  have  included  wall  shear 
stress  [lO],  Nusselt  number  Nu  and  Schmidt  number  Sc 

[ll] . None  have  been  used  to  correlate  flow  patterns  with 
corrosion  rate,  and  none  have  been  tested  for  varying  flow 
situations  . 

In  fact,  turbulence  is  not  characterized  solely  by  the 
use  of  the  Reynolds  number.  Reference  to  Figure  2 will 
quickly  and  vividly  point  out  the  possible  differences  in 
structure  that  can  readily  be  imagined  to  have  a varied 
effect  on  corrosion  including  corrosion  produce  buildup 
and  film  removal.  Both  pictures  were  obtained  by  Prandtl 
in  an  open  channel  with  aluminum  powder  sprinkled  on  the 
water  surface  [12].  In  both  pictures  the  fluid  velocity 
is  the  same;  the  only  difference  is  the  speed  of  the  camera 

li 

moving  downstream  parallel  to  the  flow.  Figure  2(a)  was 
taken  with  the  camera  moving  at  a slow  velocity.  The 
aluminum  particles  near  the  wall  appear  as  points  . In 
Figure  2(b)  , the  velocity  of  the  camera  is  about  equal  to 

that  of  the  mid-channel  flow,  making  the  wall  turbulence  ' 

appear  "streaky."  The  term  "streaky"  also  applies  to  the 
surfaces  of  the  corroded  specimens  from  the  present  experiment. 
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Figure  2.  Visualization  of  turbulent  flow  in  an  open 

channel.  (a)  Camera  speed  slow.  (b)  Camera 
speed  approximately  that  of  the  centerline 
fluid.  From  Prandtl  [12]. 
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The  implication  is  that  pure  dissolution  by  diffusion  is 
being  supplemented  by  convective  transport  of  reaction 
components  toward  and  away  from  the  metal  surface . 

The  purpose  of  the  present  investigation  was  two-fold. 
First,  it  was  a goal  to  construct  a corrosion  test  facility 
with  features  that  would  fulfill  the  requirements  for 
"smooth"  hydrodynamics,  electrochemical  monitoring  and 
visual  observation  of  the  flow.  Second,  it  was  intended  to 
experimentally  characterize  the  flow  and  pursue  possible 
correlations  with  corrosion  rate. 
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II.  DESIGN  OF  CORROSION  APPARATUS 


A.  SYSTEM  CAPABILITIES 

A photograph  of  the  corrosion  test  facility  used  in  this 
study  is  shown  in  Figure  3,  with  a schematic  shown  in  Figure 
4.  Sea  water  was  pumped  from  Monterey  Bay  into  a 19,000 
liter  (5000  gallon)  redwood  tank.  From  there,  it  was  drawn 
off  as  needed  to  fill  a 204  liter  (54  gallon)  reservoir 
drum  made  from  15  inch  PVC  pipe.  The  system  was  designed 
for  use  with  either  natural  or  artificial  sea  water,  with 
various  possibilities  for  pumping  and/or  gravity  flow. 
Specifically,  capabilities  were  provided  to: 

(1)  Fill  the  drum  from  the  tank  by  gravity. 

(2)  Fill  the  drum  with  artificial  sea  water  or  fresh  water. 

(3)  Recirculate  through  a bypass  valve. 

(4)  Recirculate  through  the  test  cell. 

(5)  Flow  through  the  test  cell  by  gravity,  directly  from 
the  tank  and  out  the  drain  line . 

(6)  Drain  the  drum  by  gravity  or  pump. 

The  redwood  tank  was  located  70  m away  from  and  5 m above  the 
test  apparatus . This  tank  was  pumped  full  through  a sand 
and  gravel  filter  1 m in  diameter  by  4 m high,  buried  in 
the  high  tidal  zone.  A Monel  base,  fitted  with  vertical 
tubes,  excluded  most  sand  and  plankton.  This  tank  was  half 
emptied  anc  refilled  once  or  twice  every  24  hours,  depending 
on  tidal  times  and  heights  and  flow  requirements . The  water 
was  filtered  again  through  a fine  mesh  paper  cartridge 
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located  between  the  tank  and  the  drum,  and  a similar  filter 
was  also  included  in  the  recirculation  loop  itself.  This 
aided  in  preventing  the  accumulation  of  foreign  matter  on 
the  anemometer  probe  that  was  used  for  turbulence  measurements. 
It  also  removed  particulates,  organisms,  corrosion  products 
and  other  substances,  such  as  colloidal  sulphur,  which  might 
be  present,  and  which  must  be  removed  in  order  that  they  do 
not  alter  the  corrosion  rate  and/or  mechanism. 

NACE  Standard  TM-01-69  gives  a preferred  minimum  volume- 
to-area  ratio  of  40  milliliters  of  solution  per  square  cen- 
timeter of  exposed  surface.  In  these  tests,  a reservoir  of 

204  liters  (54  gallons)  was  used,  giving  a ratio  of  13,600 
2 

ml/cm  . This  larger  volume: 

(1)  Reduced  the  effect  of  evaporation  or  leakage  on 
solution  properties . 

(2)  Reduced  the  effect  of  buildup  of  corrosion  products/ 
metal  ions . 

(3)  Allowed  settling  of  the  electrolyte  to  release  any 
dispersed  bubbles . 

Metal  in  the  system  was  kept  to  an  absolute  minimum,  the 
only  exception  being  the  stainless  steel  rotometer  float. 

On  tests  at  2.4  m/sec  (8  ft/sec)  with  90/10  Cu-Ni,  Ross 
found  that  additions  of  up  to  130  ppb  of  copper  ions  (maximum 
tested)  had  no  significant  effect  on  corrosion  rate  [13]. 

Cornet  has  shown  that  even  a copper  ion  concentration  of 
7 milligram  per  liter  (700  ppb)  has  no  measurable  effect 
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[14].  To  even  reach  this  level  of  contamination,  with  a 
corrosion  rate  of  6 mpy , would  require  a total  weight  loss 
of  1.43  grams  and  would  take  6130  hours  (255  days  ). 

The  NACE  standard  goes  on  to  suggest  a test  duration 
for  moderate  or  low  corrosion  rate  studies  determined  by 
the  expression: 


Test  duration  (hours) 


2000 

corrosion  rate  (mpy ) 


For  90/10  Cu-Ni,  with  a corrosion  rate  of  roughly  6 mpy, 
this  gives  14  days,  but  to  obtain  a larger  supply  of  data, 
48  hours  was  chosen  instead.  This  shorter  duration  must 
be  kept  in  mind  when  comparing  these  corrosion  rates  with 
others  in  the  literature.  As  Efird  found  with  tests  on 
90/10  at  a low  velocity  of  0.6  m/sec  (2  ft/sec),  the 
corrosion  rate  continued  to  drop  for  14  years  and  had  not 
yet  stabilized.  It  decreased  from  0.5  mpy  after  1 year  to 
0.05  after  14  years,  a ten-fold  decrease  [15]. 

B.  FLOW  LOOP 

The  flow  loop  shown  in  Figure  4 was  constructed  from 
polyvinyl  chloride  (PVC)  pipe,  fittings  and  valves.  A 
centrifugal  pump  made  of  glass  filled  epoxy  was  used,  with 
no  metal  parts  coming  in  contact  with  the  sea  water.  With 
a one  horsepower  motor,  it  provided  a maximum  flow  rate 
through  the  test  cell  of  12  gpm  at  19  pounds  force  per 
square  inch  discharge  pressure . This  corresponded  to  a 
maximum  velocity  of  7.8  m/sec  (25.6  ft/sec).  This  flow 


capability  covers  the  range  of  velocities  normally  encountered 
in  sea  water  piping  systems.  Condensors,  for  example,  in 
both  shore  based  power  plants  and  ships  at  sea,  run  from 
1 to  3 m/sec,  a compromise  between  pumping  power,  heat  trans- 
fer requirements,  and  biofouling  costs.  The  maximum  flow 
capability  of  the  present  system  corresponds  to  a speed  of 
15.6  knots,  making  higher  velocities  desirable  in  the  future. 

Flow  rate  was  controlled  by  the  pump  discharge  valve, 
a bypass  valve  and  a valve  at  the  inlet  to  the  test  cell. 

The  bypass  valve  permitted  recirculation  from  the  drum 
through  the  pump  and  back  to  the  drum.  This  allowed  for 
stabilization  of  temperature  and  oxygen  content  prior  to  the 
start  of  a test  run,  as  well  as  a bypass  for  water  dis- 
charged by  the  pump  that  did  not  flow  through  the  test  cell. 
Velocity  determinations  were  based  on  a Fisher  Porter  roto- 
meter  with  maximum  flow  capacity  of  74.2  1/min  (19.6  gal/min) . 
Calibrations  with  both  fresh  water  (specific  gravity  1.000) 
and  natural  sea  water  (specific  gravity  1.028  at  33  ppt 
salinity)  showed  an  accuracy  of  plus  or  minus  0.5  percent 
(0.06  m/sec).  Temperature  of  the  sea  water  was  periodically 
monitored  with  a thermometer  mounted  at  a 90  degree  bend 
V in  the  piping  immediately  upstream  from  the  test  cell. 

A stable  reference  electrode  was  required  from  which  all 
potentials  could  be  measured.  Saturated  calumel  has  been 
used  by  most  researchers,  but  is  quite  fragile.  The  type 
chosen  was  that  most  common  in  shipboard  cathodic  protection 
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systems  --  silver/silver  chloride.  When  a mixture  of  silver 
and  silver  chloride  interacts  with  sea  water  chloride  ions, 
a stable  reference  potential  is  generated  [16].  Unlike  the 
SCE,  the  Ag/AgCl  billet  type  electrode  could  be  easily 
mounted  in  the  piping  system  itself  as  shown  in  Figure  5 . 

The  electrode  was  screwed  into  a horizontal  "dead  leg"  at 
a 90  degree  bend  in  the  piping.  Being  horizontal,  the  well 
area  was  well  circulated  but  not  disturbed  by  turbulence 
fluctuations  or  air  pockets  that  could  give  false  readings . 

Ag/AgCl  electrodes  have  a standard  reduction  potential 
of  0.2223  volts,  referenced  to  hydrogen.  Saturated  calumel 
electrodes  (SCE)  by  comparison  are  0.2415.  Therefore,  values 
reported  in  this  paper  could  be  easily  compared  to  others  in 
the  literature  by  the  factor  0.0192  volts.  Ag/AgCl  has  the 
additional  benefit  of  relative  insensitivity  to  the  distance 
between  it  and  the  center  of  the  sample,  in  this  system  47.6 
cm  (18.75  in).  These  electrodes  can  be  periodically  re- 
juvenated with  a 1.5  volt  dry  cell  battery  and  saturated 
potassium  chloride.  A brownish  coating  of  silver  chloride 
is  deposited  on  the  silver  billet,  and  if  not  rubbed  off 
produces  consistent  results  for  several  weeks,  perhaps 
longer . 

C.  PORTABLE  TEST  CELL 

The  test  cell  itself  was  partly  modeled  after  that  used 
by  the  Francis  LaQue  Corrosion  Laboratory  in  Wrightsville 
Beach,  N.C.  and  by  Ross  [17].  Since  construction,  it  has 
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PVC  PIPE  CAP 


Figure  5.  Ag/AgCl  reference  electrode  with 
cross  section  of  holder. 


( 


also  been  found  that  Smith  used  a similar  setup  in  the  late 
sixties  [18].  The  flow  channel  had  an  internal  cross 
section  1 cm  by  1 cm  (0.394  in  square),  accurate  to  plus  or 
minus  0.0002  cm.  It  was  made  of  Plexiglas  to  prevent 
metallic  contamination  and  to  allow  viewing  of  the  flow 
during  tests.  Slip  couplings  on  each  end  permitted  easy 
removal.  As  Figures  6 and  7 show,  it  is  of  laminated  con- 
struction to  allow  disassembly  for  cleaning,  as  well  as 
substitution  of  different  Plexiglas  spacers  to  change  the 
channel  width. 

An  entrance  length  of  51  hydraulic  diameters  (four  times 
the  cro'ss  sectional  area  divided  by  the  perimeter)  was  used. 
Hinze  states  that  40  should  be  considered  a minimum  to  en- 
sure a fully  developed  velocity  profile,  with  less  entrance 
length  required  for  disturbed  entry  flow  and  high  Reynolds 
numbers,  both  of  which  are  present  here  [19].  Using  the 
empirical  formula  for  the  hydrodynamic  boundary  layer  across 
a flat  plate,  at  4 m/sec,  a layer  thickness  of  0.5  cm,  (half 
the  channel  width)  is  obtained  at  an  inlet  distance  of  33  hy- 
draulic diameters,  so  the  layers  from  opposite  walls  have 
merged  and  we  have  fully  developed  flow  over  the  specimen 
positions.  The  square  test  cell  inlet  (1.414  cm  diagonal) 
was  tapered  from  the  round  cross-section  of  the  inlet  pipe 
(1.58  cm  diameter)  to  prevent  skewing  of  the  velocity 
profile,  and  an  exit  length  of  10  hydraulic  diameters  was 
used  to  prevent  disturbance  of  this  profile  along  the 
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Figure  6.  Plexiglas  test  cell. 


Layout  of  the  Plexiglas  test  cell  showing  the 
laminated  construction  and  specimen  cavity. 


specimens.  A large  discharge  pipe  of  2.093  cm  (0.824  in) 
inside  diameter  prevented  back  pressure  disturbances,  but 
throttling  on  this  side  was  often  needed  to  prevent  cavita- 
tion across  the  sample  caused  by  the  sudden  pressure  drop. 

The  samples  were  set  into  cavities  on  one  or  both  sides 
of  the  flow  channel,  flush  with  the  Plexiglas  spacers.  Sample 
size  was  7.50  cm  by  1.40  cm  (2.95  by  0.55  in),  and  approx- 
imately 0.318  cm  (0.125  in)  thick.  Thinner  samples  would 
bow  in  the  channel  when  clamped,  while  excessively  large 
ones  would  waste  material  and  go  beyond  the  range  of  the 
available  analytical  balance  (100  grams).  It  is  considered 
that  a large  enough  exposed  area  was  used  to  reduce  the 
impact  of  any  edge  effects  on  weight  loss.  Specimens  were 
held  in  place  against  the  spacers  by  a series  of  Plexiglas 
backing  blocks  and  rubber  gaskets,  as  shown  in  Figure  8. 

Electrical  connection  to  the  specimens,  when  needed,  was 
provided  by  brass  screws  tapered  to  a dull  point  and 
threaded  into  the  backing  blocks . Semicircular  cavities 
along  the  side  of  one  specimen  cavity  allowed  easy  removal 
of  even  a tight  specimen  with  a special  pair  of  pliers. 

The  specimen  on  the  other  side  could  then  by  pushed  out, 
ensuring  no  metal  removal  that  might  invalidate  weight 
loss  measurements . 

Flow  through  the  test  cell  was  turbulent  (Reynolds 

i 

number  based  on  diameter  was  greater  than  2300).  For  the 

three  velocities  used  in  the  present  experiment,  the  , 

corresponding  Reynolds  numbers  were: 
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Figure  8.  Arrangement  for  mounting  specimens  in  test 
showing  parts  used  for  polarization  plots  . 
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In  an  attempt  to  change  the  character  of  the  turbulence, 
i.e.  eddy  size,  energy  distribution  and  turbulence  intensity 
at  a given  average  velocity,  a variable-mesh  grid  network 
arrangement,  shown  in  Figure  9,  was  placed  upstream  of  the 
specimen  positions  in  the  flow.  The  mesh  was  made  with 
0.015  cm  ( 0 . 006  in)diameter  three  pound  test  fishing  line 
strung  in  a Plexiglas  frame,  with  four  different  size 
meshes  . 

D.  CHOICE  OF  ELECTROLYTE 

1 . Synthetic  Versus  Natural  Sea  Water 

Both  synthetic  (artificial)  and  natural  (fresh)  have 
been  used  in  corrosion  studies.  On  the  one  hand,  artificial 
sea  water  is  reprodacibly  formulated  [20],  with  pH  and  con- 
ductivity constant.  When  synthetic  sea  water  is  used, 
correlation  of  results  is  often  simplified.  Studies  using 
synthetic  sea  water  can  be  particularly  effective  in  inves- 
tigating differences  between  various  metals  or  the  effect  of 
a single  isolated  parameter  on  corrosion.  On  the  other  hand 
some  effects  that  take  place  in  a natural  environment  are 
tremendously  altered  or  obliterated  when  using  synthetic 
sea  water  in  a laboratory  situation.  The  use  of  natural  sea 


Figure  9 


Arrangement  of  Plexiglas  mesh-strung  block 
placed  in  entry  section  of  test  cell. 
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water  are  often  several  times  greater  than  in  synthetic 
sea  water.  Admittedly,  natural  sea  water  is  an  extremely 
complex  mixture.  While  it  varies  in  content  and  properties 
around  the  world,  well  circulated  supplies  that  are  well 
aerated  and  free  of  pollution  can  be  expected  to  produce 
similar  results  in  different  locations. 

2 . Supply  of  Natural  Sea  Water 

a.  Description  of  Monterey  Bay,  California 

As  shown  in  Figure  10,  Monterey  Bay  is  an  open 
embeyment  located  on  the  central  California  coast.  Rich 
in  flora  and  fauna,  the  mouth  of  the  bay  is  37  km  (20 
nautical  miles)  wide.  Bisecting  the  bay  almost  symmetrically 
is  the  Monterey  Submarine  Canyon,  one  of  the  world's  largest. 
Its  head  comes  within  0.3  km  (0.16  miles)  of  the  shoreline, 
with  depths  of  200  m (656  ft)  only  one  mile  from  shore,  and 
increases  in  depth  seaward  to  865  m (2840  ft),  over  half 
a mile. 

Good  circulation  of  the  bay  is  assured  by  the 
current  system.  The  California  Current  flows  slowly  south- 
ward, transporting  low  temperature,  low  salinity,  highly 
oxygenated,  nutrient  rich  Subarctic  water.  The  California 
Countercurrent  flows  nearer  the  surface,  proceeding  north- 
ward, transporting  warm,  high  salinity,  oxygen  poor,  nutrient 
rich  Equatorial  Pacific  water  [21  ]. 
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Figure  10.  Map  of  Monterey  Bay  showing  location  of  sea 
water  supply  to  corrosion  test  facility. 
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From  February  to  September,  strong  winds  push  the 
surface  water  offshore,  allowing  subsurface  water  to  replace 
it;  this  is  a process  known  as  upwelling.  This  makes 
coastal  water  low  in  temperature  and  oxygen  content  and 
high  in  salinity  and  nutrient  content.  During  the  oceanic 
period,  from  September  to  November,  the  dense,  previously 
upwelled  water  sinks,  and  offshore  surface  water  flows 
on  shore.  Surface  temperatures  reach  an  annual  maximum, 
salinity  levels  off,  and  nutrient  concentrations  rise. 

To  complete  the  cycle,  during  the  Davidson  Current  period, 
a local  surfacing  of  the  California  Countercurrent , from 
November  to  February,  water  from  the  south  converges  against 
the  coast  and  sinks.  Surface  temperature  decreases  abruptly, 
but  subsurface  temperature  increases.  Salinity  is  low  due 
to  surface  runoff  and  rain  (in  a normal  season). 

General  circulation  within  the  bay  appears  to  be  that 
of  onshore  movement  in  the  central  region  and  counter- 
clockwise and  clockwise  flow  in  the  northern  and  southern 
areas  respectively.  The  residence  time  of  water  between 
the  central  bay  and  the  northern  and  southern  ends  varies 
from  3 to  8 days,  depending  on  time  of  year.  The  coastline 
is  moderately  polluted,  subject  to  agricutural  as  well  as 
industrial  runoff,  but  with  good  circulation  [22]. 

Sea  water  for  this  experiment  was  obtained  through 
a sand  and  gravel  filter  buried  in  the  high  tidal  zone  at 
a location  36°  36.5'  north  latitude  and  121°  53.0'  west 


longitude.  The  nearest  extensive  hydrographic  data  [23] 
is  available  from  a sampling  site  approximately  0.1  miles 
away,  where  the  water  depth  is  16  m (52  ft).  For  the  period 
from  January  1972  to  April  1973,  Figure  11  shows  the  varia- 
tion in  surface  temperature.  Figure  12  shows  the  salinity, 
which  was  computed  from  the  conductivity  [24].  Figure  13 
shows  the  observed  dissolved  oxygen  concentration  (measured 
by  the  Winkler  method  of  titration) . All  three  figures  are 
for  samples  taken  from  the  surface  of  the  bay.  The  limited 
variation  of  all  three  variables  is  significant.  In  particular, 
an  annual  variation  in  temperature  of  only  5°C  makes  it  pos- 
sible to  perform  experiments  year  around,  whereas  those 
done  in  some  other  parts  of  the  world  might  yield  results 
strongly  dependent  upon,  but  inseparable  from,  the  tempera- 
ture change.  As  examples,  the  Francis  LaQue  Corrosion 
Laboratory  has  reported  seasonal  variations  in  temperature 
and  salinity  of  40-81°F  and  32-37  ppt  [25],  while  the  Ocean 
City  Research  Corporation  has  reported  variations  of  32-80°F 
and  29-36  ppt  [8]  . 

3 . Sea  Water  Variables 
a.  Marine  Fouling 

In  terms  of  corrosion,  the  differences  between 
artificial  and  natural  sea  water  may  become  more  pronounced 
if  marine  organisms  influence  results . During  low  velocity 
tests,  some  metals  are  prone  to  fouling  by  marine  life, 
including  slime  film,  barnacle  cyprids  (which  cannot  attach 
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Figure  11.  Variation  in  surface  water  temperature  in  Monterey  Bay  (36°  36.6'  north 
latitude,  121  53.1'  west  longitude).  Data  from  Broenkow  [23]. 


Monterey  Bay  (36  36.6'  north 

from  Broenkow  [23]. 


at  velocities  greater  than  0.5  knots)  and  tubeworms  to  name 
the  most  prevelant  [26].  In  stagnant  or  polluted  water, 
aerobic  forms  can  lower  the  oxygen  content,  allowing  sulfate 
reducing  bacteria  to  take  control,  producing  hydrogen  sul- 
fide gas  and  dark  colored  water  [18].  A 1969  National  Bureau 
of  Standards  Report  on  "Microbial  Corrosion"  cites  an  increase 
in  corrosion  rate  after  one  month  of  six  times  that  before 
inoculation  with  a particular  strain  [27].  Rapid  reproduction 
of  certain  species  of  bacteria  has  been  reported  by  Bott, 
where  a single  slime  cell  3 ^m  in  diameter  increased  by  an 
order  of  magnitude  in  size  in  only  4 hours. 

Increasing  the  fluid  velocity  makes  slime  attachment 
more  difficult.  The  slime  layer  becomes  thinner,  denser  and 
more  tightly  adhered  to  the  s’-rface  [28].  Figures  14  and  15 
show  the  decrease  in  fouling  rate  (increase  in  fouling 
resistance)  with  increased  velocity  for  both  90/10  Cu-Ni 
and  titanium ;"the.  Cu-Ni  was  more  resistant  to  fouling  due 
to  its  toxicity  [29].  Below  1.8  m/sec,  the  fouling  of  90/10 
was  very  sporadic,  periodically  sluffing  off,  while  at 
higher  velocities  it  fouled  at  a fairly  steady  rate. 

Long  term  tests  may  similarly  be  influenced  by  this 
fouling.  Photographs  of  90/10  Cu-Ni  by  Efird  after  14  years 
showed  only  light  fouling,  while  other  "non-toxic"  metals 
may  become  encrusted  with  several  pounds  per  square  foot  in 
only  one  year  [15].  Fouling  is  much  more  variable  from 
place  to  place  in  the  oceans  than  is  corrosion. 
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Figure  14.  Fouling  rates  for  titanium  and  90/10  Cu-Ni 
(CDA  706).  From  Ritter  [29], 
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Figure  15 . Effect  of  velocity  on  fouling  rate 
Ritter  [29 ] . 
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b.  Dissolved  Oxygen 

Oxygen  is  known  to  be  a major  controlling  vari- 
able in  corrosion.  As  it  was  not  controlled  in  this  experi- 
ment, its  variability  must  be  discussed.  Two  sources  supply 
oxygen  to  the  sea  water.  First  is  the  sea/air  interface, 
affected  by  wind,  waves,  tide  and  time  of  the  year  (upwelling, 
for  example) , where  water  becomes  saturated  to  varying 
degrees.  Generally,  cold  water  has  higher  oxygen  concentra- 
tions than  warm  water  because  oxygen  saturation  values  at 
a given  atmospheric  partial  pressure  depend  inversely  on 
temperature  and  salinity,  with  salinity  being  the  dominant 
factor  [22].  As  an  example  of  the  variation  possible,  Scott 
found  surface  values  of  dissolved  oxygen  within  the  bay  it- 
self during  the  upwelling  period  to  vary  from  6.3  to  8.6 
ppt  from  one  location  to  another  (approximately  7.0  nearest 
our  source  of  water)  [30]. 

A second  supply  of  oxygen  is  biochemical.  During 
photosynthesis,  phytoplankton  remove  nutrients  from  the  water 
and  release  oxygen,  while  during  decomposition  nutrients 
are  released  and  oxygen  consumed  [2l].  This  biological  pro- 
duction is  controlled  by  solar  radiation,  upwelling,  tem- 
perature and  depth  of  the  mixed  layer,  and  nutrient  supply 
[22]. 

Horne  has  reported  that  where  the  concentration 
of  nutrients  is  high,  such  as  phosphates  and  nitrates,  the 
oxygen  content  tends  to  be  low  [31].  The  three  great  oceans 
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of  the  world  tend  to  have  characteristic  dissolved  oxygen 
profiles.  The  Atlantic  receives  an  influx  of  new,  cold, 
oxygen  rich  water  from  both  polar  regions.  In  contrast, 
the  Indian  Ocean  and  the  northern  boundaries  of  the  Pacific 
are  landlocked,  their  waters  being  older  and  relatively 
depleted  in  oxygen.  There  the  circulation  is  poorer  and 
more  organic  material  has  been  oxidized  in  them  for  a longer 
period  of  time. 

One  additional  generalization  about  the  oxygen  con- 
tent of  the  seas  is  its  variability  with  depth.  The  surface 
layer  is  well  mixed  with  oxygen  concentration  near  saturation 
and  uniform  down  to  the  thermocline . Below  this  layer,  a 
subsurface  maximum  may  occur  as  a result  of  oxygen  produc- 
tion by  photosynthesis.  Below  this,  in  situ  processes  such 
as  oxidation  of  organic  matter  tend  to  reduce  the  oxygen 
content.  Yet  deeper,  it  again  increases  due  to  the  influx 
of  young,  dense,  cold,  oxygen  rich  water  that  has  sunk 
down  from  the  polar  regions  [31]. 

This  discussion  serves  to  point  out  not  only  the 
variability  possible  between  different  waters  of  the  world, 
but  even  between  water  in  the  same  location  during  different 
seasons  or  different  years. 

The  saturation  solubility  of  oxygen  in  sea  water  can 
be  calculated  from  the  equation  given  by  Truesdale,  a cubic 
function  of  temperature  and  a linear  function  of  salinity. 
However,  even  easier  is  the  use  of  his  reduced  formula: 
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C ' = a - bS 
s 

where  ' is  the  oxygen  solubility  in  ppm,  S is  the  salinity 
in  ppt , and  a and  b are  found  from  his  tables  using  the 
temperature  (zero  to  40°C)  , at  0.1°  intervals  [32].  The 
trend  is  for  a decrease  in  solubility  with  increasing  salt 
concentration.  The  dissolved  oxygen  content  in  this  experi- 
ment was  determined  using  a Hach  Model  0X-2P  Dissolved 
Oxygen  Test  Kit.  It  is  based  on  the  Winkler  titration 
method  and  accurate  to  0.2  ppm.  Values  encountered  ranged 
from  5.2  to  6.8  ppm,  while  the  corresponding  saturation 
levels  were  7.9  and  6.6  ppm. 

The  effect  of  oxygen  on  corrosion  rate  is  many- 
faceted.  It  can  act  to  depolarize  cathodic  areas,  oxidize 
metal  ions  in  the  electrolyte  to  a higher  valence  or  form 
a protective  film.  Since  corrosion  is  an  electrochemical 
process,  it  can  be  divided  into  an  oxidation  and  reduction 
reaction.  The  rate  of  corrosion  is  determined  by  the 
"slower  step."  For  copper  alloys  in  sea  water  (alkaline), 
the  reduction  reaction  is: 

O2  + 2H2O  + 4e" ►40H" 

Reinhart  found  only  a slight  increase  in  corrosion  rate 
as  oxygen  was  increased  from  0.4  to  5.75  ppm  [33].  Recently 
MacDonald  and  Syrett  [34]  have  extended  this  range  up  to 
26.3  ppm;  although  this  concentration  is  more  than  double 
the  value  that  would  be  found  in  natural  waters,  the  dramatic 
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increase  in  corrosion  rate  highlights  not  only  the  importance 
of  oxygen  itself,  but  also  of  any  other  variable  that  in- 
creases the  quantity  of  dissolved  oxygen  that  reaches  the 
metal  surface.  MacDonald  and  Syrett  hypothesize  that  the 
increase  in  corrosion  rate  at  high  oxygen  content  is  due  to 
a positive  shift  of  the  corrosion  potential  beyond  a critical 
value,  where  a different  reaction  mechanism  becomes  possible, 
greatly  enhancing  film  breakdown  for  certain  copper  alloys 
[35].  This  occurs  despite  the  microscopic  observation  that 
a thicker  film  is  formed  in  the  more  highly  oxygenated  water. 
They  found  that  a high  oxygen  levels  the  reaction  rate  is 
under  anodic  control,  where  the  corrosion  rate  is  determined 
by  ionic  or  electronic  transport  through  the  surface  film.  At 
very  low  oxygen  contents,  cathodic  control  would  dominate, 
with  the  corrosion  rate  determined  by  the  transport  of 
oxygen  to  the  surface. 

c.  Temperature,  Salinity,  Conductivity,  pH 

Temperature  also  has  an  effect  on  corrosion. 

While  in  general,  corrosion  kinetics  tend  to  increase  with 
temperature,  this  effect  is  counteracted  by  the  decreasing 
solubility  of  oxygen  in  water  at  higher  temperatures.  Also 
as  temperature  increases,  the  diffusion  coefficient  of  oxygen 
increases,  which  tends  to  increase  corrosion  rate,  but  this 
change  is  tempered  by  the  concurrent  increase  in  the  thick- 
ness of  the  diffusion  boundary  layer  which  acts  to  decrease 

corrosion  rate.  Also,  an  increase  in  temperature  can 

; 

increase  the  rate  of  marine  fouling. 
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Temperatures  within  the  laboratory  itself  occassion- 
ally  ranged  from  8 to  20°C  in  a 12  hour  period.  However, 
ordinarily  only  a 5°C  variation  was  experienced,  which 
led  to  a fluctuation  of  up  to  3°C  in  the  flow  loop  sea  water 
temperature.  The  original  sea  water  temperature  was  also 
raised  by  the  heating  effect  of  the  pump,  reaching  a steady 
state  condition  of  6 to  8°C  above  room  temperature.  Recir- 
culating the  water  for  6 to  8 hours  prior  to  a test  run 
reduced  the  temperature  rise  during  a run  from  10  to  only 
2°C.  Peterson  and  Gehring  have  recently  tested  a new  flow 
channel  where  the  temperature  variation  was  even  more  sig- 
nificant, 21.6  to  41.0°C  [36]. 

Salinity,  expressed  in  ppt , is  the  percentage  of 
solid  materT al  present,  and  is  roughly  a measure  of  the 
total  salt  c./ntent.  It  is  usually  determined  by  measuring 
the  electricax  conductivity  of  the  sample,  and  then  con- 
verting with  appropriate  tables  from  Cox  [24].  The  specific 
conductance  of  sea  water  varies  almost  linearly  with  tempera- 
ture and  electrolyte  concentration,  i.e.  salinity.  At  a 
salinity  of  33  to  34  ppt,  it  is  0.035,  0.040  and  0.045  mho/cm 
for  10,  15  and  20°C  respectively  [3l].  Conductivity  in  this 
experiment  ranged  from  0.0308  to  0.0371  mho/cm. 

In  natural  sea  water,  pH  is  controlled  primarily  by 
the  amount  of  dissolved  CO2  in  the  water  [31].  Photosynthetic 
activity  and  solar  radiation  are  controlling.  For  surface 
water,  equilibrium  with  the  atmosphere  usually  assures  a 
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high  pH.  Normal  values  for  sea  water  are  7.5  to  8.4,  slightly 
alkaline.  Interactions  of  CO^,  ammonia  and  hydrogen  sulfide 
(formed  on  the  death  and  subsequent  decay  of  fouling  organ- 
isms) can  also  produce  local  perturbations  in  pH.  It  can 
thus  vary  seasonally  and  daily.  During  this  study  the 
range  was  from  7.68  to  7.92,  measured  with  a 4403  Electro- 
mark Digital  Analyzer  accurate  to  0.01.  Porte  has  shown 
that  even  at  temperatures  as  high  as  40°C,  the  corrosion 
rate  is  independent  of  pH  for  the  range  5 to  8 , increasing 
slightly  for  pH  less  than  4 and  decreasing  slightly  for  a 
pH  above  10  [9  ] . 
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III.  EXPERIMENTAL  PROCEDURE 


A.  POTENTIODYNAMIC  POLARIZATION  PLOTS 

As  a metal  corrodes,  it  assumes  a particular  potential. 
This  potential  can  be  measured  directly,  but  the  corrosion 
current  cannot.  The  corrosion  current  is  equal  to  the 
sum  total  of  all  the  current  that  is  flowing  between  the 
local  anodes  and  cathodes  on  the  metal  surface.  It  can  be 
measured  indirectly  using  a potentiostat , by  controlling 
the  potential  of  the  corroding  metal  and  measuring  the  cell 
current.  The  resultant  curve  is  a polarization  plot.  The 
potentiostat  can  drive  the  potential  from  the  least  noble 
potential  desired  through  the  equilibrium  corrosion  potential 
Ecorr  most  noble  potential  desired.  The  upper  and 

lower  portions  of  the  curve  correspond  to  the  anodic  and 
cathodic  reactions  respectively.  The  slopes  of  the  linear 
portions  of  these  two  arms  are  the  so-called  Tafel  slopes. 

The  intersection  point  of  the  two  Tafel  lines  gives  the  equil- 
ibrium corrosion  potential  (E  ) and  the  corrosion  current 
(Icorr^  for  that  particular  period  in  time.  The  plot  is 
made  at  a relatively  slow  scan  rate  of  1 mV/sec,  so  the 

corrosion  current  found  is  really  an  average  over  the  15 
•» 

minutes  required  to  run  the  plot. 

B.  HOT  FILM  AIJEMOMETRY 

In  order  to  characterize  the  flow  situation  within  the 
test  cell,  measurements  were  made  with  a hot  film  anemometer 
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mounted  as  shown  in  Figure  16.  The  detecting  element  of  the 
probe  is  a very  fine  metal  film  on  a glass  substrate  0.025  mm 
(0.001  in)  in  diameter  and  0.5  mm  (0.020  in)  long.  The  film 
is  heated  by  an  electric  current,  but  the  flowing  water  tends 
to  cool  the  film,  lowering  its  electrical  resistance. 
Therefore,  the  current  needed  to  maintain  the  film  at  a con- 
stant temperature  can  be  related  to  the  fluid  velocity  (in 
a linear  manner  with  proper  circuitry) . The  film  itself 
is  coated  with  a thin  layer  of  quartz  to  protect  it  from 
abrasion  by  suspended  particles,  and  the  joints  and  supports 
of  the  probe  are  insulated.  Figure  17  shows  an  unused 
probe  on  the  right  and  a probe  used  for  less  than  one  hour 
on  the  left.  While  the  fouling  was  easily  removed  by  ultra- 
sonic cleaning,  it  did  alter  the  calibration,  tended  to  give 
suspicious  results  and  led  to  premature  failure.  Provision 
of  an  additional  filter,  in  the  flow  loop  itself,  solved 
the  fouling  problem,  which  had  existed  even  when  fresh  tap 
water  was  used. 

In  this  study,  a miniature  probe  was  used.  The  small 
frontal  area  of  the  probe  shaft,  1.5  mm  (0.06  in)  minimized 
the  increased  local  velocity  that  results  when  any  protru- 
berence  reduces  the  local  flow  channel  cross  sectional  area. 
Hinze  has  shown  that  a small  sensing  wire  is  needed  for  low 
velocities  and/or  small  eddy  measurements  [19].  However, 
for  measurements  close  to  the  wall,  even  this  small  obstacle 
to  flow  is  significant.  Fluid  flows  between  it  and  the  wall. 
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Figure  16.  Plexiglas  block  (fits  into  specimen  cavity)  with 

anemometer  probe  and  cross  section  of  probe  holder 


Figure  17,  Cylindrical  hot  film  probes.  Probe  on 

left  used  less  than  one  hour  in  unfiltered 
fresh  water.  Probe  on  right  unused  (22x). 
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accelerating  and  separating.  Only  with  an  external  device, 
such  as  a laser,  can  undisturbed  measurements  be  made. 

When  flow  characterization  measurements  were  being  made, 
a Plexiglas  block  was  fitted  into  the  specimen  cavity  on 
one  side,  making  three  sides  of  the  flow  channel  completely 
smooth  for  their  entire  length.  The  hot  film  probe  was 
mounted  in  the  other  specimen  cavity  through  a hole  in  the 
block.  It  was  mounted  near  the  center  to  avoid  any  dis- 
turbances caused  by  the  leading  edge  of  the  block  or  by  the 
corners  of  the  channel,  and  hopefully  would  give  data  that 
represented  the  typical  flow  structure  over  the  corroding 
metal  surface . 

C.  CORROSION  TESTING 

1 . Sample  Preparation  and  Handling 

The  material  used  as  samples  in  this  study  was  90/10 
copper-nickel  (CDA  Alloy  706)  hot  rolled  5/8  inch  plate  with 
soft  temper.  Chemical  composition  is  given  in 

Table  2,  with  the  mechanical  properties  given  in  Table  3. 
Figure  18  shows  the  microstructure  of  a polished  specimen, 
etched  with  a mixture  of  1 part  hydrogen  peroxide  (30  per 
cent),  19  parts  distilled  water,  and  10  parts  ammonium  hy- 
droxide (specific  gravity  0.90)  Figure  18(a)  shows  the 
average  grain  size,  which  varied  with  distance  from  the 
rolled  surface.  Figure  18(b)  shows  certain  crystallographic 
planes  within  the  grains . It  is  apparent  that  a great  deal 
of  distortion  still  exists  within  the  metal,  a condition 
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TABLE  2 


MIL-C-15726  E Amendment  2 CDA  706 


Element 


Composition  (per  cent) 


Copper 

87.40 

Nickel 

10.40 

Iron 

1.50 

Manganese 

0.49 

Zinc 

0.13 

Tin 

0.02 

Silicon 

0.02 

Lead 

0.01 

Phosphorous 

0.013 

Sulphur 

0.001 

Chemical  composition  of  90/10  Cu-Ni . 


TABLE  3 


Yield  Strength 
Tensile  Strength 
Elongation 
Hardness 
Dens ity 


Mechanical  properties 


.1068  MN/m^  (24.0  ksi) 
. 2358  MN/m"^  (53.0  ksi) 
40  per  cent 
Rockwell  B-25 
8.93  g/cm^ 


f 90/10  Cu-Ni. 
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Microstructure  of  annealed  90/10  Cu-Ni, 
polished  and  etched.  (a)  Average  grain 
size.  (lOOX)  (b)  Crystallographic  planes 
within  grains.  (lOOX) 


that  would  increase  the  corrosion  rate.  The  horizontal 
shadows  in  both  photographs  are  termed  structural  streaks , 
revealed  by  the  etchant  as  a result  of  structural  heter- 
ogeneities within  the  plate. 

Samples  for  polarization  plots  done  in  a beaker  were 
0.953  cm  (0.375  in)  high  by  0.953  cm  (0.375  in  ) diameter 
right  circular  cylinders  with  an  exposed  surface  area  of 
6.467  cm^  (1.002  in^).  A 3-48  hole  2 mm  deep  was  used  to 
attach  the  sample  to  the  electrode  holder.  Samples  for  the 
flow  channel  were  machined  so  a transverse  face  was  exposed. 
Samples  were  polished  with  3/0  grade  emery  paper.  The  first 
four  sets  of  samples  were  run  with  polishing  lines  parallel 
to  the  flow  to  minimize  erosion  by  accentuated  turbulence, 
and  get  accurate  weight  loss  data.  Subsequent  samples  were 
run  with  polishing  lines  perpendicular  to  the  flow,  to 
allow  differentiation  between  streamwise  corrosion  features 
that  result  from  preferential  dissolution  of  the  lines 
(stress  concentration  sites,  or  oxygen  deficiency  at  their 
bases),  and  those  features  caused  by  the  flow  itself.  The 
samples  were  then  wiped  clean  with  Kim-Wipes,  degreased 
with  ethyl  alcohol,  and  air  dried.  The  dried  samples  to  be 
used  for  weight  loss  determinations  were  weighed  on  a 
Sartorius-Werke  A6  balance  to  an  accuracy  of  0.1  milligrams. 

Generally,  two  samples  were  run  simultaneously  on 
opposite  sides  of  the  flow  channel.  One  was  used  for  weight 
loss  and  the  other  for  SEM  study.  The  SEM  sample  was  cut 
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into  three  pieces  2.5  by  1.4  cm  (0.98  by  0.55  in),  any  one 
of  which  would  fit  easily  onto  a large  SEM  stub  of  diameter 
3.15  cm  (1.24  in).  Usually,  only  the  cut  section  contain- 
ing the  leading  edge  was  used,  as  it  contained  the  general 
morphology  present  in  the  other  sections.  This  also  allowed 
observation  of  crevice  corrosion  or  leading  edge  flow  effects, 
if  present. 

After  the  tests,  specimens  for  weight  loss  were 
cleaned  by  dipping  for  five  minutes  in  a mixture  of  one 
part  HCl  (specific  gravity  1.18),  5 parts  sulfuric  acid 
(specific  gravity  1.84)  and  4 parts  distilled  water,  as 
recommended  in  the  ASTM  Standard  Gl-72(37).  This  was 
followed  by  light  scrubbing  in  distilled  water  with  a fine 
bristle  toothbrush,  rinsing  in  ethyl  alcohol  and  air  drying. 
The  effect  of  this  cleaning  solution  in  removing  meral  was 
tested  on  a channel  specimen  and  found  to  be  negligible 
within  the  accuracy  of  the  scale  used.  After  exposure, 
the  specimens  were  weighed  to  the  same  accuracy  of  plus  or 
minus  one  mg.  The  corrosion  rate  was  calculated  from  the 
weight  loss  by  the  formula: 


Corrosion  rate  = 


constant  X weight  loss 


area  X time  X metal  density 


2 

with  weight  loss  in  grams,  area  in  cm  , time  in  hours  and 

3 

density  in  g/cm  . To  obtain  the  corrosion  rate  in  mpy , 
the  constant  is  3.45  X 10^,  while  8.76  X 10^  will  yield 
the  value  in  Mmpy . 
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2 . Polarization  Resistance  Measurements 

The  rate  at  which  a metal  corrodes  can  vary  with 
time;  often  the  variation  is  significant.  Monitoring  this 
change  can  give  useful  information  concerning  which  corrosion 
mechanisms  are  operating  and  the  structure  of  the  protective 
film.  The  instantaneous  corrosion  rate  and  corrosion  current 
can  be  found  from  the  polarization  resistance,  i.e.  the  re- 
sistance of  the  metal  surface  to  a change  in  applied  potential 
and  equals : 

R - E (applied) 
p ” A i (measured ) 

where  A i is  the  change  in  current  per  unit  surface  area  of 
the  specimen.  The  potential  of  the  specimen  (working 
electrode)  as  compared  to  the  reference  electrode  is  varied 
by  passing  a suitable  current  between  the  working  and  the 
counter  electrode. 

The  corrosion  current  is  related  to  R by  the 

P 

Stern-Geary  formula  [38]:  ' 

_ ^a  fic 

^corr  " 2 . 303  ( j3a  + j3cT 

The  sensitivity  of  this  formula  to  the  Tafel  slopes  can  be 
seen  from  the  factor  ^a  ^c/(  fia.  + /Sc).  For  fia  - - 

0.12  Volts/decade,  the  factor  is  0.06,  while  for  ^a  - ^c  = 
0.60,  the  factor  is  0 . 3,  which  is  higher  by  a factor  of  5. 

The  corrosion  rate  in  mpy  is  then  calculated  from  Faraday's 
Law  by  the  formula: 


„ , o Equivalent  weight  . 
“Py  = DensTty  ^ ^ 


corr 


Macdonald,  Syrett  and  Wing  have  obtained  good  corre- 
lation between  polarization  resistance  measurements  and 
weight  loss  by  calculating  the  equivalent  weight  using  a 
composition-averaged  atomic  weight  (M)  and  change  in  oxida- 
tion state  (n)  [34]: 


M = M^  + X^.M... 
Cu  Cu  Ni  Ni 


The  mole  fractions  and  Xj^^^  are  calculated  from  the 

chemical  composition  data  given  in  weight  per  cent.  Oxida- 
tion states  for  the  copper  and  nickel  are  assumed  to  be  +1 
and  +2  respectively,  although  others  may  also  be  operative. 

In  the  present  experiments,  the  test  cell  was  instru- 
mented to  conduct  these  measurements  during  several  corrosion 
runs.  A Princeton  Applied  Research  Model  173  Potentiostat 
with  a Model  175  Programmer  supplying  the  required  input 
and  monitoring  was  used,  with  a Hewlett  Packard  Model  7040A  X-Y 
Recorder  for  output.  To  disturb  the  metal  surface  as  little 
as  possible,  an  imposed  voltage  of  only  ^ 12  mV  was  applied. 
Because  this  "tickling"  may  affect  the  subsequent  corrosion 
rate  as  well  as  measure  the  instantaneous  rate,  the  possibility 
of  it  affecting  the  weight  loss  measurements  was  avoided  by 
using  one  specimen  for  polarization  resistance  (and  SEM 
analysis  in  the  as-corroded  condition) , and  a separate  one 
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for  weight  loss  (and  SEM  analysis  of  the  clean  surface).  ' 

Measurements  were  made  at  varying  intervals  throughout  the 

run,  with  more  measurements  in  the  early  stages,  when 

corrosion  rate  was  higher.  A relatively  slow  scan  rate 

(0.5  mV/sec)  was  used  to  assure  accurate  results. 

Macdonald  has  recently  shown  that  a sizable  error  in 
corrosion  rate  determination  can  result  from  even  a scan 

rate  only  ten  times  this  value,  as  well  as  a sizable  hy-  ; 

steresis  loop  [34].  The  scan  was  done  first  in  the  anodic 

direction  and  then  the  cathodic  direction,  forming  a counter-  ; 

clockwise  hysteresis  loop.  It  should  be  noted  that  these 
plots  were  made  with  the  potential  on  the  ordinate  and  the 
current  on  the  abscissa;  with  the  axes  reversed,  the  loop 
would  be  clockwise.  The  maximum  hysteresis  observed  was 
3.5  /jlA  after  3 minutes.  By  10  minutes,  it  dropped  to  0.5 
fxA  and  after  15  minutes  was  not  detectable.  Using  a 
forward  and  reverse  sweep  also  gave  two  curves  that  could 
be  averaged  to  get  a better  value  of  A I . As  Walter  points 
out,  this  curve  will  be  nearly  linear  only  if  the  Tafel 
slopes  are  equal  [39].  This  was  the  case  in  the  present 
study . 

To  increase  the  accuracy  of  results,  the  effect 
of  IR  drop  must  also  be  considered.  This  accounts  for  the 
resistance  of  the  electrolyte  between  the  reference  and 
working  electrodes.  As  natural  sea  water  has  quite  high 
conductivity,  IR  drop  was  neglected  in  this  study.  However, 
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it  was  found  that  doubling  the  exposed  surface  area  (by 
connecting  the  potentiostat  to  both  channel  specimens 
simultaneously)  did  not  result  in  a doubling  of  AI . This 
may  have  been  due  to  neglecting  an  appreciable  IR  drop, 
in  which  case  calculated  corrosion  rates  would  be  too  low. 

One  additional  aspect  of  polarization  resistance 
measurements  that  has  thus  far  lacked  attention  is  the  deter- 
mination of  exposed  surface  area.  In  calculating  mpy,  the 
original  "flat"  surface  area  is  customarily  used.  However, 
most  surfaces  actually  have  more  metal  exposed  than  this  would 
indicate,  due  to  surface  flaws,  undulations  and  machining 
marks.  Even  a 3/0  grade  emery  finish  gives  an  rms  roughness 
of  approximately  10  microinches.  The  surface  roughness 
often  becomes  exaggerated  as  corrosion  proceeds  . For  a 
metal  surface  that  is  50  per  cent  covered  with  hemispherical 
pits  2 ptm  in  diameter,  the  surface  area  exposed  to  the 
corrosive  liquid  is  actually  17  per  cent  greater  than  a 
perfectly  flat  surface.  The  corrosion  rate  calculated, 
by  either  polarization  resistance  measurements  or  weight 
loss,  using  the  smaller  area  may  then  be  too  high  by  perhaps 
half  that  value:  8.5  per  cent.  If  the  corrosion  is  in  the 
form  of  hemispherical  troughs,  as  in  the  present  study,  the 
increase  in  surface  area  rises  to  29  per  cent. 
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IV.  RESULTS  AND  DISCUSSION 


A.  LEADING  EDGE  ATTACK 

Some  investigators,  such  as  Syrett  [10]  and  Ross  [40], 
have  reported  preferential  dissolution  of  the  leading  edge 
of  specimens,  and  attributed  this  to  the  establishment  of 
the  diffusion  mass  transfer  (or  oxygen  concentration) 
boundary  layer.  Cornet  has  also  mentioned  this  entrance 
effect,  where  the  so-called  oxygen  tension  and  the  copper 
ion  concentration  are  changing  [11].  To  check  this  affect, 
split  samples  were  placed  end-to-end  and  run  at  three  dif- 
erent  velocities  in  the  flow  channel.  These  specimens 
were  made  from  full-size  90/10  Cu-Ni  coupons  cut  in  half, 
making  a total  of  4 metal  specimens,  two  on  each  side  of 
the  flow  channel.  Thus,  for  2 of  the  4 specimens,  the 
leading  edge  was  preceded  by  Plexiglas,  while  the  other  2 
further  downstream  were  preceded  by  metal.  Results  of  these 
tests  are  shown  in  Table  4.  In  all  cases  except  one,  the 
trailing  section  actually  had  a greater  weight  loss  than 
the  leading  section,  differing  by  a maximum  of  only  5.2 
per  cent.  The  corrosion  rate  was  consistently  higher  for 
specimens  from  one  side  of  the  flow  channel  with  a maximum 
of  18.5  per  cent  difference.  The  slight  difference  in  exposed 
area  of  the  specimens  on  the  opposing  sides  was  taken  into 
account  by  measuring  the  specimen  length  and  the  channel 
width  on  both  sides . The  reason  then  must  be  a flow  phenom- 
enon, though  the  nearest  upstream  pipe  bend  is  1 m away,  and 
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TABLE  4 


Velocity 

Channel 

Corrosion 

Leading 

Rate  (mpy) 
Trailing 

Difference 

(m/ sec ) 

Side 

Specimen 

Specimen 

(mpy ) 

(per  ce 

4,0 

1 

6.40 

6.74 

0 . 34 

6.2 

4.0 

2 

5 . 33 

5.48 

0.15 

2.8 

6 . 0 

1 

14.96 

14.99 

0.03 

0.2 

6.0 

2 

13.61 

13.35 

0 .26 

1.9 

7.8 

1 

18 . 31 

18.51 

0.20 

1.1 

7.8 

2 

16.79 

16.82 

0 .03 

0.2 

Split  specimen  direct  weight  loss  data 
taken  at  various  fluid  velocities  (48 
hour  tests ) . 
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lA 


I 


f;  the  nearest  nonuniformity  inside  the  flow  channel  is  51  cm 

away.  Microscopic  examination  showed  no  preferential  attack 
p whatsoever  at  either  leading  edge  in  48  hours. 

I Three  possible  explanations  exist  besides  the  diffusion 

boundary  layer  concept  for  why  other  investigators  may  have 
found  this  preferential  attack.  First  is  the  problem  of 
uneven  mating  surfaces.  The  experiments  of  Locke  have 
shown  that  even  for  a metal  tightly  bonded  in  an  epoxy  mold 
and  subsequently  sanded  smooth,  the  metal  bexng  harder  stays 
slightly  raised  above  the  epoxy  [41].  This  results  in  pre- 
ferential attack  at  the  interface,  due  to  the  establishment 
of  the  metal  corner  as  'a  local  anode . His  tests  were  done 
in  beakers  with  stagnant  artificial  sea  water.  In  a flow 
situation,  the  effect  of  a metal  surface  sticking  out  in 
the  flow  path  would  be  obviously  amplified  by  the  induced 
turbulence  and  erosion  "tearing"  at  the  leading  edge. 

The  opposite  case  of  an  indented  or  recessed  metal 
surface  could  easily  occur  in  an  arrangement  such  as  that 
described  in  this  paper.  A particle  of  dirt  between  the 
' sample  and  the  channel  spacers,  or  an  exposed  surface 

slightly  rounded  at  the  edge  by  sanding  on  a soft  table 
would  result  in  a stream  of  cavitation  bubbles  emanating 
from  the  leading  edge.  This  was  noticed  in  trial  runs, 
but  avoided  in  all  reported  runs  by  careful  sample  polishing, 
cleanliness  in  the  sample  cavity,  and  external  clamping  of 
the  specimens  against  the  channel  spacerSi. 


67 
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The  third  difficulty  may  involve  a crevice  at  the 
leading  edge,  yielding  cavitation,  increased  turbulence 

and  crevice  corrosion  (deficiency  of  oxygen  resulting  in  j 

accelerated  corrosion) . This  was  avoided  by  machining  the 

specimens  slightly  oversize  and  then  hand  sanding  with  3/0  < 

grade  emery  paper  to  fit  tightly.  Microscopic  examination  | 

1 

revealed  a square  leading  edge  as  opposed  to  the  slightly  j 

rounded  edge  caused  by  the  milling  operation.  All  of  these  j 

defects  would  go  unnoticed  in  an  arrangement  where  viewing 

of  the  flow  was  not  possible  either  because  the  sample  i 

itself  was  moving  or  because  the  apparatus  was  not  transparent. 

After  removal  of  the  specimen,  an  explanation  of  the  cause 

of  any  preferential  corrosion  would  be  highly  speculative.  ' 

These  results  challenge  the  design  notion  that  insertion 
of  an  entry  length  of  test  material  is  required.  It  may, 
in  fact,  damage  the  experiment,  if  it  involves  recirculation 
of  a fixed  volume  of  electrolyte,  by  multiplying  the  metal 
ion  concentration.  Additionally,  a corrodable  entry 
length  may  change  the  flow  conditions  with  time  as  pitting 

or  surface  roughness  developes  . j 

Further,  in  regard  to  the  mass  transfer  boundary  layer,  i 

1 

the  T-7C-5  Task  Group  of  the  National  Association  of  Corro- 
sion Engineers  (NACE) , in  a 1977  review  of  high  velocity  sea  ' 

water  corrosion,  used  the  following  equation  for  the  boundary  ’ 

layer  thickness  when  the  hydrodynamic  boundary  layer  also 

I 

starts  at  the  same  location  [42]:  j 
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However,  this  shows  that  it  does  not  stabilize  rapidly  at 

4/5 

all,  but  increases  in  proportion  to  X .A  "dummy  section" 

ahead  of  the  test  section,  would  still  result  in  S in- 

m 

creasing  across  the  specimen.  As  Leumer  has  pointed  out 
from  a review  of  the  literature  on  mass  transfer,  use  of 
this  formula  in  turbulent  flow  is  inappropriate  [43].  In 
laminar  flow,  while  the  formula  applies,  it  shows  that  use 
of  an  entry  length  would  not  prevent  a concentration  gra- 
dient across  the  specimen;  the  layer  thickness  does  not 
stabilize  . 


B.  POLARIZATION  BEHAVIOR 

Figures  19  and  20  show  the  standard  potentiodynamic 
polarization  curves  obtained  for  CDA  Alloy  706  (right  cir- 
cular cylinder  specimen  in  a standard  corrosion  cell  flask) 
in  synthetic  and  natural  sea  water  respectively.  Comparison 
of  the  two  plots  shows: 
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(1)  The  equilibrium  potential  is  0.005  volts  more 
negative  (cathodic)  in  natural  sea  water  due  to  the  larger 
cathodic  Tafel  slope. 

(2)  The  corrosion  current  (after  ten  minutes  immersion) 
is  36.5  p,A  in  synthetic  sea  water,  54  /xA  in  natural  sea 
water,  a 39  per  cent  difference. 

Figure  21  shows  a plot  that  was  also  done  in  natural 
sea  water,  with  stirring  from  a magnetic  agitator.  However, 
this  plot  was  performed  in  the  portable  Plexiglas  flow 
channel  cell,  isolated  from  the  rest  of  the  system  by  rubber 
stoppers.  This  test  arrangement  was  considerably  different: 

(1)  Rectangular  sample  with  an  exposed  surface  area  of 

2 . 2 . 

7.5  cm  (1.16  in  ) vice  a right  circular  cylinder  with  an 
area  of  6.47  cm^  (1.00  in^). 

(2)  Ag/AgCl  reference  electrode  located  47.6  cm  (18.75  in) 
from  the  sample  with  agitation  only  near  the  sample  itself, 
vice  5.0  cm  (2.0  in)  away  with  agitation  around  all  electrodes. 

(3)  Platinum  (counter)  electrode  was  positioned  1.0  cm 
from  the  sample  (working  electrode),  vice  a carbon  electrode 
5.0  cm  (2.0  in)  away. 

The  plots  are  shifted  slightly  along  the  current  axis 
due  to  the  different  exposed  areas  of  the  two  types  of 
samples.  Comparison  of  the  corrosion  current  density  ( yuA/cm  ) 
shows  them  to  be  quite  similar  — 8.34  in  the  standard  cell 
versus  8.40  in  the  channel,  using  a cathodic  Tafel  slope  of 
0.4  volts  per  decade.  The  equilibrium  potentials  of  speci- 
mens in  the  flask  and  flow  channel  are  -0.23  and  -0.24  volts 
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respectively,  and  the  overall  shapes  of  the  curves  are  like- 
wise very  similar.  Thus,  this  flow  channel  electrode  arrange- 
ment has  been  proven  to  give  accurate  results,  in  terms 
of  it  reproducing  the  data  in  a carefully  controlled  lab- 
oratory test  following  ASTM  Standard  G5-72  [44]. 

Figure  22  is  another  polarization  plot  performed  in  the 
flow  channel  approximately  one  month  later.  The  shapes  of 
the  two  curves  are  quite  different,  although  only  three 
known  variables  had  changed; 

(1)  A carbon  counter  electrode  positioned  39.4  cm  (15.5  in) 
downstream  from  the  specimen  was  used  in  place  of  the  platinum 
electrode  positioned  directly  across  the  channel.  However, 
similar  plots  were  also  run  at  this  time  with  the  platinum 
electrode  and  the  same  anomaly  observed. 

(2)  The  sea  water  temperature  was  16.5  vice  14.5°C. 

(3)  The  test  was  performed  with  a flow  velocity  provided 
by  the  pump  of  6 m/sec,  vice  the  low  speed  gravity  flow. 

However,  tests  at  the  same  period  with  even  an  isolated  test 
cell  produced  curves  similar  to  Figure  22.  Peterson  and  Gehring 
have  recently  experienced  a similar  change  with  tests  of  low 
carbon  steel  [36].  As  in  their  situation,  the  discrepancy 
cannot  be  explained  by  a variation  in  oxygen  content. 

Plots  as  in  Figure  22  were  obtained  for  varying  velocities 
with  the  following  shift  toward  more  negative  (more  noble) 
corrosion  potentials  (after  10  minutes  immersion): 
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ion  plot  obtained  in  flow  channel  for  90/10 
natural  sea  water  at  14.5°C,  gravity  flow  of 


Velocity  (m/sec) 


corr 


(Volts  vs.  Ag/AgCl) 


2.0  -0.265  J 

1 

4.0  -0.270  < 

6.0  -0.273 

7.8  -0.276  j 

i 

Additionally,  the  slope  of  the  anodic  portion  of  the  curve 

increased  slightly  with  increased  velocity,  showing  an 

increase  in  the  corrosion  current  (I  ).  The  switch- 

corr 

back  shown  in  Figures  19,  20,  and  21,  at  a potential  of 
about  0.15  volts,  corresponds  to  a change  in  corrosion 
mechanism:  the  formation  of  black,  film,  probably  cuprous 
oxide  (CuO).  This  oxide  layer  was  observed  while  making 
all  3 of  these  plots,  but  not  during  construction  of 
Figure  22. 


C.  EFFECT  OF  VELOCITY  ON  CORROSION  RATE 

Velocity  effects  vary  from  alloy  to  alloy  depending  mostly 
on  their  respective  film  formation  characteristics.  While 
some  alloys  form  very  tenacious  protective  films,  others 
form  easily  perturbed  film  and  still  others  none  at  all. 

Copper  alloys  typically  form  films  that  withstand  attack 
remarkably  well  at  low  velocities;  a slight  increase  in 
corrosion  rate  with  velocity  is  usually  seen,  due  to  enhanced 
oxygen  transport.  However,  as  velocity  increases,  a point 
is  eventually  reached  where  the  film  breaks  down  and  typically 
begins  to  be  locally  removed.  Danek  lists  this  so-called 
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breakaway  velocity  as  4.5  m/sec  (15  ft/sec),  for  90/10 
Cu-Ni  [5], while  Syrett  showed  it  to  also  be  a function  of 
the  test  system  geometry  [lO].  It  is  undoubtedly  a function 
of  both  average  velocity  and  geometry,  but  more  fundamentally 
a function  of  the  flow  structure  very  near  the  metal  surface. 
At  very  high  velocities,  another  regime  may  be  encountered 
where  a film  cannot  form  and  corrosion  rate  stabilizes,  but 
this  is  well  above  the  range  of  interest  of  the  present  study 
and  of  present  day  technological  applications. 

Due  to  lack  of  a fundamental  model,  corrosion  rates  in 
this  study  are  simply  plotted  versus  average  velocity,  as  is 
typical  in  "velocity-effects"  studies.  In  comparing  corrosion 
rate  results  from  various  sources,  the  means  of  estimating 
the  average  velocity  which  is  cited  must  be  kept  in  mind. 
Usually,  in  tube  flow  experiments  the  average  flow  velocity 
is  calculated  from  the  quantity  of  water  moving  through  the 
tube  divided  by  the  cross  sectional  area  of  the  tube;  this 
gives  a velocity  that  is  less  than  the  centerline  velocity, 
although  the  difference  decreases  with  increased  flow  rate. 

On  the  other  hand,  systems  that  employ  a rotating  or  revolving 
electrode  use  the  relative  free  stream  velocity  which,  stir- 
ring effects  neglected,  gives  a velocity  higher  than  the 
average . 

The  corrosion  rates  at  varying  velocities  from  the  present 
study  are  shown  in  Table  5 and  Figure  23.  Direct  weight  loss 
measurements  were  converted  to  mpy  for  all  runs.  For 
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Corrosion  rates  of  90/10  Cu-Ni  as  calculated  directly  from 
weight  loss  data  and  from  polarisation  resistance  data  (48 
hour  tests ) . 


comparison  of  calculated  corrosion  rates,  polarization  re- 
sistance data  was  also  taken  for  2 of  the  runs,  and  converted 
to  mpy . By  integrating  the  area  under  the  curve  of  polariza- 
tion resistance  versus  time,  an  average  corrosion  rate  was 
calculated,  for  comparison  with  the  direct  weight  loss 
determinations.  In  these  two  methods,  the  following  errors 
and  uncertainties  are  present: 

(1)  Flow  adjustment,  as  read  from  the  rotometer. 

(2)  Varying  temperature,  salinity,  pH,  dissolved  oxygen 
content . 

(3)  Measurement  of  exposed  surface  area  of  the  specimens. 

(4)  Electrical  and  hydrodynamic  noise  in  the  polarization 
resistance  plots. 

(5)  Weighing  of  specimens  during  use  of  the  analytical 
balance  or  due  to  incomplete  removal  of  corrosion  products. 

(6)  Correction  for  IR  drop  was  not  included. 

The  calculations  of  corrosion  rate  (mpy)  from  polariza- 
tion resistance  measurements  (Table  5)  were  based  on  Tafel 
slopes,  /8^  and  , of  0.12  volts  per  decade,  considered  an 
average  for  many  metal/electrolyte  systems.  Due  to  the 
shapes  of  the  polarization  curves  in  Figures  21  and  22, 
actual  Tafel  slopes  would  be  impossible  to  measure  accurately. 
The  curve  in  Figure  21  shows  slopes  of  approximately  0.08, 
although  for  the  cathodic  portion  this  estimation  is  ex- 
tremely arbitrary.  Use  of  these  values  would  yield  calculated 
corrosion  rates  33  per  cent  lower  than  those  found  using  0.12. 
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At  the  other  extreme.  Figure  22  shows  no  linear  portion;  a 
slope  of  0.6  (or  even  higher)  could  be  selected.  Using  0.6 
would  increase  the  calculated  corrosion  rates  400  per  cent 
above  the  0.12  t-ised  values.  Even  0.2  would  yield  a 67 
per  cent  increase.  Thus,  unless  accurate  Tafel  slopes  are 
known,  polarization  resistance  data  is  most  useful  in  observ- 
ing trends,  i.e.  relative  changes  in  corrosion  rate  (with 
time  and/or  velocity)  , as  opposed  to  measuring  absolute 
values . 

All  of  the  corrosion  rates  calculated  in  Table  5 from 
direct  weight  loss  measurements  are  higher  than  those  found 
by  other  investigators  (Figure  1)  especially  at  6 and  7.8 
m/sec.  Since  tests  in  the  present  study  were  only  of  48 
hour  duration  (vice  times  of  up  to  months  in  other  studies), 
and  since  the  corrosion  rate  decreases  with  time,  these  higher 
initial  values  are  to  be  expected. 

These  corrosion  rates  have  been  compared  to  those  found 
by  Leumer  using  a circling  foil  apparatus  and  synthetic 
sea  water  [43].  Natural  sea  water  is  a more  aggressive  en- 
vironment, as  shown  by  the  higher  value  of  I for  natural 

sea  water  (Figures  19  and  20).  In  spite  of  this  fact,  the 
corrosion  rates  in  the  flow  channel  averaged  approximately 
20  per  cent  lower  than  those  for  the  circling  foil  at  4 m/sec. 
This  suggests  that  some  additional  factor  of  corrosion  was 
experienced  on  the  circling  foil,  such  as  due  to  erosion  by 
entrained  bubbles,  or  due  to  higher  turbulence  intensity. 


At  6 m/sec,  the  situation  was  reversed,  with  the  flow 
channel  giving  corrosion  rates  roughly  50  per  cent  higher 
than  the  foil.  While  the  explanation  for  this  reversal 
is  not  known,  it  does  support  the  conviction  that  better 
modeling  of  the  corrosion  process  is  required. 

Only  one  of  the  possible  inlet  meshes  was  used  in  the 
present  study,  with  the  mesh  oriented  at  a 45  degree  angle 
to  the  channel  walls,  and  spacing  between  filaments  of  2.36 
mm.  Reference  to  Table  5 shows  that,  at  least  for  a fairly 
low  velocity  (4  m/sec),  insertion  of  a mesh  at  the  test 
cell  entrance  had  no  effect  on  the  corrosion  rate  (as 
calculated  by  weight  loss).  Either  the  change  in  turbulence 
near  the  metal  surface  had  no  effect  on  corrosion,  or  the 
change  in  turbulent  structure  introduced  by  the  mesh  dis- 
appeared by  the  time  the  flow  reached  the  sample ; the  latter 
is  believed  to  be  the  case. 

Figure  24  indicates  the  dramatic  effect  of  temperature 
on  corrosion  rate.  After  the  first  few  hours,  as  the  pro- 
tective film  develops,  the  corrosion  rate  would  be  expected 
to  decrease.  However,  in  the  case  shown  in  Figure  24,  the 
*:emperature  rise  dominates  the  corrosion  rate,  causing  it  to 
rease  over  60  per  cent.  That  the  two  variables  follow  so 
- . ; . certainly  interesting,  but  the  data  is  too  limited 
' . = ive.  Syrett  has  summarized  the  literature  on  the 

[lO].  Various  investigators  have  found 
• ’ - * . increase,  decrease  or  be  unaffected 
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Figure  24.  Variation  of  temperature  anci  corrosion  rate  during  test 
at  6 m/sec.  Sample  area  = 7,5  cm^ , = 0.020  Volts. 


by  an  increase  in  temperature.  Syrett  attributes  this 
variation  to  varying  dissolved  oxygen  solubility  and 
differences  in  the  uniformity  and  tenacity  of  the  corro- 
sion product  film.  For  the  case  presented  here,  the 
oxygen  content  only  dropped  from  6.8  to  5.6  ppm,  with  the 
initial  drop  to  5.8  in  the  first  24  hours.  The  inter- 
dependence is  obviously  a complex  one. 

D.  SURFACE  MORPHOLOGY 

Figures  25(a)  and  (b)  are  SEM  photomicrographs  taken  at 
llOOX  and  2300X  respectively  of  the  cleaned  surface  of  a 
specimen  which  had  been  corroded  for  48  hours  at  7.8  m/sec 
(25.6  ft/sec).  Comparison  with  Figure  18(b)  will  point 
out  the  preferential  dissolution  along  certain  crystallo- 
graphic planes  within  the  base  metal  grains . The  pre- 
ferentially attacked  areas  are  those  which  first  become 
covered  with  corrosion  products.  As  Tait  points  out, 
oxygen  concentration  cells  then  form  and  subsequently  ini- 
tiate pits  [45].  Pores,  ruptures  and  flaws  in  the  film 
ensure  uneven  diffusion  of  oxygen,  allowing  oxygen  concentra- 
tion cells  to  keep  these  new  pits  growing. 

Metallic  copper  that  is  removed  from  the  surface  goes 
either  into  the  corrosion  product  film  (Cu^O),  or  into  solu- 
tion as  ionic  copper  (toxic  to  some  marine  life),  or  forms 
a precipitate  in  the  solution.  In  tests  with  90/10  Cu-Ni 
in  0.1  per  cent  NaCl,  Popplewell  found  that  for  the  first 
100  hours,  more  than  80  per  cent  of  the  metal  lost  was  used 
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Figure  25.  SEM  photomicrographs  of  cleaned  surface  of 
specimen  exposed  for  48  hours  at  7.8  m/sec. 
Direction  of  fluid  flow  is  horizontal,  pol- 
ishing lines  vertical,  (a)  llOOX  (b)  2300X 
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to  form  a protective  film.  Only  after  1000  hours  had  the 
film  thickness  stabilized,  with  all  the  metal  subsequently 
lost  going  into  solution  [46].  Samples  in  the  present 
experiment,  run  at  4 and  6 m/sec  (13.1  and  19.7  ft/sec), 
showed  no  visible  signs  of  corrosion  product  buildup  after 
48  hours  time;  the  corroded  surface  was  essentially 
featureless  with  small,  random  groups  of  shallow  pits. 

Weighing  the  specimen  before  and  after  post-exposure  clean- 
ing showed  only  a negligible  weight  change  (less  than  0.01 
mg/cm  ).  However,  the  samples  run  at  7.8  m/sec  formed  a 
"streaky"  corrosion  product  layer  on  top  of  trough-like  pits, 
as  shown  in  Figure  26.  Of  the  weight  loss  during  corrosion, 
8.9  to  10.8  per  cent  was  found  to  be  in  the  corrosion  product 
that  adhered  to  the  metal  surface.  While  wet,  the  corrosion 
product  was  rust-colored  and  easily  wiped  off.  When  dry, 
it  was  tightly  adhered,  purple  to  the  naked  eye  and  green 
under  the  microscope.  In  order  to  get  a view  of  the  corro- 
sion pits  at  a further  stage  of  development,  but  without 
running  a longer  test,  a specimen  was  coupled  to  a platinum 
electrode  of  the  same  exposed  surface  area  located  directly 
across  the  channel.  The  resultant  uncleaned  surface  of  the 
specimen  exposed  in  this  galvanic  couple  is  shown  in  Figure  27 
The  pits  are  larger  than  in  an  uncoupled  specimen,  as  expected 
but  there  is  a lack  of  corrosion  product  buildup. 

The  streaks  of  corrosion  products  and  pitting  parallel 
to  the  direction  of  flow  indicate  a strong  dependence  on 
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Figure  26.  As-exposed  surface  of  specimen  exposed  for  48 

hours  at  7.8  m/sec.  Flow  direction  is  horizon- 
tal, polishing  lines  vertical.  (200X) 


Figure  27.  As-exposed  surface  of  specimen  exposed  for  48 
hours  at  7.8  m/sec,  coupled  to  a platinum 
electrode  on  the  opposite  side  of  the  channel. 
Flow  direction  is  horizontal,  polishing  lines 
vertical.  (lOOX) 
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flow  phenomena . These  streaks  are  spaced  more  closely  than 
the  structural  streaks  referred  to  earlier  (Figure  18(a)) 
and  they  are  more  periodic  in  distribution  and  more  uniform 
in  size.  It  is  doubtful  that  they  are  connected  directly 
to  the  turbulence  concept  of  streaks  of  fluid  lifting  off 
the  wall  (the  process  of  bursting  described  later  in  this 
paper)  as  this  bursting  is  a more  random  phenomenon  than 
the  distribution  of  corrosion  attack  seems  to  exhibit.  The 
periodicity  and  elongation  in  the  line  of  flow  do  appear 
to  be  strongly  influenced  by  the  flow,  perhaps  in  a manner 
similar  to  the  production  of  period  sand  ripples  on  the  ocean 
floor  (aiso  discussed  in  a later  section). 

The  green  color  of  the  corrosion  product  indicates  it  to 
be  a copper  hydroxychloride  [34].  However,  potential-pH 
diagrams  which  show  the  thermodynamically  stable  forms  of 
the  metal  and  its  corrosion  products,  show  that  conditions 
in  the  present  experiment  fall  within  the  stability  range  of 
cuprous  oxide.  No  corrosion  product  was  observed  to  collect 
on  the  in-line  filter. 

E.  STRUCTURE  OF  TURBULENCE 

1 . General  Concepts 

For  the  practicing  engineer,  the  more  of  the  details 
of  turbulence  that  can  be  neglected  the  better.  However, 
in  this  section  it  is  necessary  to  go  into  some  detail  to 
give  the  reader  at  least  a "gut  feeling"  of  which  turbulence 
parameters  may  be  important  to  corrosion  and  why.  We  will 
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specifically  restrict  ourselves  to  a discussion  of  the 
turbulent  (vice  laminar)  flow  regime.  As  Prandtl  pointed 
out  in  1934,  many,  if  not  most,  problems  of  interest  in  real 
life  fall  into  this  category  [l2].  For  example: 

(1)  Flow  through  a smooth  pipe  with  inside  diameter  2.54 

cm  (1.00  in)  will  be  turbulent  greater  than  2300)  for 

any  flow  rate  over  0.65  liters/min  (0.17  gal/min) , i.e. 
average  velocity  above  0.11  m/sec  (0.35  ft/sec). 

(2)  Flow  through  a 15.24  cm  (6.0  in)  smooth  pipe  will  be 
turbulent  above  0.11  liters/min  (0.03  gal/min),  0.02  m/sec 
(0.06  ft/sec ) . 

(3)  For  a given  velocity  of  1 m/sec,  the  maximum  pipe 
diameter  that  would  maintain  laminar  flow  would  be  0.24  cm 
(0.09  in),  giving  a maxim.um  flow  rate  of  only  0.27  liters/min 
(0.7  gal/min) . 

(4)  Flow  along  a ship  from  a point  1 m from  the  bow  to 
the  stern  would  be  turbulent  (Re^  greater  than  5 X 10^)  for 
all  velocities  greater  than  0.51  m/sec  (only  1.01  knots). 

(5)  For  a given  velocity  of  5.08  m/sec  (10  knots),  the 
laminar  portion  of  the  flow  along  the  ship  would  be  confined 
to  the  first  51  cm  (20  in). 

Figure  28  shows  the  two  flow  regimes  as  well  as  flow  that 
cannot  clearly  classified  as  either  (so-called  transition 
flow)  . 

In  corrosion,  we  are  concerned  with  wall  (vice  free) 
turbulence,  i.e.  turbulence  generated  and  continuously 
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affected  by  friction  forces  at  the  wall,  as  opposed  to  that 
produced  by  interacting  flow  of  fluid  layers  with  different 
velocities  in  the  absence  of  walls.  In  wall  turbulence, 
the  small  scale  structure  is  much  less  uniform  than  it  is 
in  free  turbulence.  The  turbulent  boundary  layer  consists 
of  flow  in  a continuous  state  of  transition.  A few  general- 
izations about  wall  turbulence  include: 

(1)  It  is  a random  phenomenon,  but  with  a quasi- 
permanancy  and  quasi-periodicity  in  time  and  space  (i.e. 
repeats  in  time  and  space,  but  not  periodically  at  one 
place  in  time  nor  at  one  time  in  space). 

(2)  It  is  characterized  by  a strong  diffusive  nature. 

(3)  Dissipation  of  turbulence  occurs  due  to  viscous 
forces,  mainly  in  the  region  of  the  smallest  eddies,  with 
kinetic  energy  converted  to  heat. 

(4)  It  is  anisotropic,  i.e.  the  mean  velocity  shows  a 
gradient  associated  with  the  shear  stress  . 

(5)  It  can  be  generated  only  by  velocity  gradients  in 
the  flow. 

(6)  Streamlines  lose  their  identity. 

(7)  The  laminar  sublayer  becomes  very  thin  and  decreases 
with  increased  velocity. 

(8)  The  average  velocity  approaches  the  free  stream  or 
centerline  velocity,  with  the  mean  velocity  gradient  changing 
rapidly  near  the  wall.  This  is  shown  in  Figure  29  at  inter- 
mediate Reynolds  numbers.  The  profile  becomes  even  more 
"square"  with  increasing  Reynolds  numbers. 
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29.  Velocity  profile  across  an  open  channel 
with  Re,  = 5600  (o)  and  8200  (+).  From 
Eckelmann  [48]. 
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(9)  Turbulence  measurements  are  often  normalized  with 

wall  parameters . The  friction  or  shear  velocity  u^  is 

1/2 

defined  by  ( '’’w/p)  , and  is  computed  from  the  slope  of 

the  mean  velocity  profile  near  the  wall.  Distance  from  the 
wall  is  then  non-dimensionalized  as  y+  = yu^/v  and  velocity 
normalized  as  u"*”  = uu^/p  with  dimension  of  sec”^ . 

For  flat  plate  or  pipe  flow,  the  velocity  profile  near 
the  wall  is  given  by  the  universal  velocity  distribution  law: 

+ 1 , + . 
u = £ In  y + c 

Very  near  the  wall,  in  the  viscous  region,  this  reduces  to 
u"*”  = y^ . Plotted  on  semi-log  paper.  Figure  30  shows  this 
relationship.  Varying  constants  have  been  found  by  investigators. 
For  example,  Wallace  [49]  found  1/k  = 2.5,  while  Eckelmann 
[48]  got  2.65.  C is  a function  of  both  surface  roughness 
and  Reynolds  number;  typical  values  obtained  by  the  above 
investigators  were  5.5  and  5.9. 

Figure  31(a)  shows  an  instantaneous  trace  of  the  axial 
flow  velocity.  Coupled  with  Figures  2 and  28,  these  indicate 
somewhat  the  complexity  of  turbulence.  On  a macroscale, 
turbulence  is  a disordered  phenomenon.  When  repeating 
velocity  patterns  are  viewed  at  a fixed  point,  they  are  never 
sufficiently  periodic  to  be  clearly  identifiable.  This 
irregularity  makes  it  impossible  to  describe  all  the  details 
of  the  motion  in  terms  of  space  and  time  coordinates.  The 
approach  then  is  to  describe  it  by  laws  of  probability,  in- 
dicating mean,  average  or  rms  values. 
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Figure  30.  Turbulent  velocity  distribution 

at  Re,  = 7150.  From  Wallace  [49]. 
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Figure  31.  Burst  identification.  (a)  Axial  velocity  u^ . 

(b)  Signal  from  spectrum  analvzer.  (c)  Time 
derivative  du,/dt.  (d)  Spatial  derivative 
du^ / dr.  From  Heidrick,  Banerjee  and  Azad  [50] 
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Burst/Sweep  Phenomenon 

The  overall  turbulence  models  that  have  been  proposed 
usually  consist  of  two  parts: 

(1)  An  organized  component  with  a quasi-periodic 
motion  of  larger  eddies,  perhaps  modeled  by  wave  theory. 

(2)  A more  random  type  of  background  turbulence, 
modeled  as  small  eddies  or  noise. 

The  precise  mechanism  which  triggers  wall  turbulence  is  not 
fully  understood.  What  is  known  is  the  sequence  of  events 
that  occur  in  the  turbulent  region.  The  scenario  based  on 
visual  studies  near  the  wall  by  Kim,  Kline  and  Reynolds 
reads  somewhat  like  a movie  script  [51].  They  observed  that 
"streaks"  of  slow-moving  and  fast-moving  fluid  side  by  side 
flowed  along  the  wall-,  randomly,  a low  speed  streak  lifted 
away  from  the  wall,  developing  into  vortices  which  migrated 
from  the  wall,  and  then  broke  up  into  more  chaotic  motion; 
these  periods  of  intense  activity  interspersed  with  relatively 
quiescent  periods  have  been  termed  "bursts"  and  "sweeps." 

For  a surface  phenomenon  like  corrosion,  this  model  of 
streaks  lifting  away  from  the  wall,  disintegrating  and  smash- 
ing back  into  the  wall  has  obvious  implications . Where  the 
corrosion  rate  is  controlled  by  the  buildup  of  a protective 
film,  a more  rapid  replenishment  of  oxygen  convected  directly 
to  the  surface  would  build  ip  the  film  faster  and  perhaps 
make  it  thicker.  Either  event,  the  movement  of  fluid  away 
from  the  surface  or  its  return  could  interupt  the  protective 
film,  exposing  fresh  metal. 
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Offen  views  this  phenomenon  of  low  velocity  fluid 
erupting  into  the  high  velocity  region  as  akin  to  separation 
caused  by  a temporary  local  adverse  pressure  gradient  [52]. 
This  pressure  oscillation  would  be  caused  by  the  arrival  of 
a sweep  (circulatory  motion)  as  it  passes  over  a newly 
forming  streak.  This  correlates  with  Wallace's  observations 
that  sweeps  precede  bursts,  followed  by  new  sweeps  further 
downstream  [49].  The  fluid  returning  to  the  wall  spreads 
out  sideways  and  is  quickly  retarded  by  the  strong  viscous 
forces  near  the  wall.  Thus  the  streak-lifting  seems  to  be 
triggered  by  large  scale  disturbances  already  present  in  the 
flow;  as  Naraynanam  has  pointed  out,  this  is  not  entirely 
a wall  phenomenon  at  all  [53]. 

Alternatively,  Michalke  explained  bursting  as  a result 
of  a locally  unstable  laminar  shear  layer  which  had  lifted 
up  and  rolled  into  two  counter-rotating  vortices  that  roll 
around  each  other  and  coalesce  [54],  Figure  32  shows  a 
computer  analysis  by  Acton  of  such  a mechanism  [55].  In 
this  pairing,  a vortex  from  a previous  burst  upstream  would 
pass  over  a vortex  associated  with  a wall  streak,  with 
portions  of  it  moving  toward  the  wall  to  provide  fluid  for 
the  next  streak.  As  these  vortices  would  not  usually  be 
vertically  aligned,  complex  motions  would  result. 

Figure  31(a)  shows  how  the  axial  velocity  u^  may  vary 
with  time.  Only  by  passing  the  unprocessed  signal  through 
a spectrum  analyzer.  Figure  31(b),  or  by  taking  the  time 
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Figure  32 


Development  of  a vortex  sheet  at 
varying  instants  in  time . From 
Acton  [55]. 
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derivative  of  the  signal  before  filtering,  Figure  31(c), 
or  by  using  the  filtered  spatial  derivative  of  the  velocity. 
Figure  31(d),  do  bursts  become  apparent. 

Bark  has  hypothesized  that  the  large  scale  eddies  are 
excited  by  the  small  scale  ones  through  the  bursting  mechanism. 
Periodically,  these  large  eddies  would  create  conditions 
suitable  for  initiation  of  new  bursts  [56].  Landahl  conversely 
proposed  that  bursting  was  due  to  the  focusing  and  trapping 
of  small  secondary  instability  waves  on  the  large  scale 
eddies  [57]  . 

Whatever  the  mechanism,  Bremhorst  has  shown  that  the 
erupted  low  momentum  fluid  diffuses  through  much  of  the 
boundary  layer,  so  that  this  is  not  at  all  a totally  localized 
event  [58].  Sabot  observed  the  bursts  to  originate  across 
most  of  the  boundary  layer,  and  then  travel  across  the  flow 
approximately  equal  to  their  streamwise  extent  [59];  he 
found  that  the  ejected  vortex  filaments  may  persist  across 
the  entire  boundary  layer,  this  being  an  important  basic 
mechanism  responsible  for  turbulent  transport  controlling 
heat,  mass  and  momentum  transfer.  Sabot  proposed  that  the 
large  scale  bursting  is  also  responsible  for  transport  of 
shear  stress,  with  the  fine  scale  structure  possibly  imbedded 
in  the  bursts,  forming  the  high  frequency  part. 

A burst  then  is  characterized  by  the  outward  ejection  of 
low  speed  fluid  in  the  form  of  small  jets  with  high  frequency 
energy,  where  the  fluctuating  streanwise  velocity  u^  is  low 
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or  decreasing.  A sweep  exhibits  the  inward  rush  of  fluid 
toward  the  wall,  when  uj^  is  high  and  increasing.  Decelera- 
tion seems  to  promote  turbulence,  while  acceleration  has 
a damping  effect.  Most  of  the  production  of  turbulence 
occurs  during  these  bursts.  However,  in  a fully  turbulent 
region,  they  may  be  obscured  by  background  turbulence. 

Increased  Reynolds  number  corresponds  to  increased  burst 
rate,  and  therefore  increased  energy  production.  The  corro- 
tion  rate  also,  in  general,  increases  with  velocity,  perhaps 
closely  tied  to  the  increased  burst  rate.  This  burst  rate 
would  determine  how  often  parcels  of  fresh  oxygen  were  de- 
livered right  to  the  doorstep  of  the  corroding  local  anode . 

A similar  mechanism  has  also  been  observed  in  fully 
developed  pipe  flow  by  Corino  and  Brodkey  [60],  Their  analysis 
pointed  to  several  stages : 

(1)  Local  deceleration  of  a large  portion  of  the  flow 
near  the  wall . 

(2)  This  slower  fluid  forced  to  accelerate  slightly 
above  the  wall  by  faster  moving  fluid. 

(3)  The  resultant  strong  shear  layer  near  the  wall  leading 
to  a small  parcel  of  low  velocity  fluid  being  rapidly  ejected 
outward . 

(4)  Followed  by  a sweeping  motion  of  high  speed  fluid 
from  upstream  toward  the  wall. 
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Together,  these  events  (bursts  and  sweeps)  produce  a 
stress  considerably  greater  than  the  total  Reynolds  stress 
(shear  stress  due  to  inertial  forces).  Wallace  has  shown 
localized  regions  where  dissipation  is  several  orders  of 
magnitude  greater  than  the  average  at  the  wall  [49].  Eckelmann 
has  found  instantaneous  values  of  Reynolds  stress  with  peaks 
30  times  the  mean  value  in  the  viscous  sublayer  itself  [48]  . 
Hinze  found  instantaneous  values  of  shear  stress  60  times  the 
local  value  just  prior  to  and  during  a burst  [l9]  . The  vortices 
can  contribute  to  the  shear  stress  only  when  they  suffer 
significant  deformation  in  cross  section  or  lateral  movement, 
which  occurs  most  strongly  during  coalescence  [6l]  . As  the 
width  of  the  eddies  increases,  so  does  the  shear  stress  [52]. 

If  wall  shear  stress  is  to  be  considered  to  have  an 
effect  on  corrosion,  it  is  these  fluctuation  peaks  that  are 
most  important.  Using  pressure  drop  measurements,  Macdonald 
[34]  has  calculated  the  average  wall  shear  stress  at  1.62 
m/sec  (5.3  ft/sec)  to  be  0.00154  pounds  force  per  square  inch. 
Using  the  Moody  diagram  for  pipe  roughness,  at  6 m/sec  in  the 
square  flow  channel  of  the  present  study,  the  calculated 
values  for  range  from  0.014  for  smooth  pipe  to  0.048  for 
a relative  roughness  (peaks  divided  by  pipe  diameter)  of  .05. 
With  fluctuations  of  two  orders  of  magnitude  in  the  shear 
stress,  it  would  indeed  seem  possible  that  strongly  affects 
corrosion  product  buildup  and  even  film  removal.  Shear  stress 
considerations  are  different  in  pipe  flow  and  external  flow, 
since  the  velocity  gradients  perpendicular  to  the  wall  vary. 
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The  wall  region  in  pipe  flow  is  often  called  the  constant 


stress  layer,  due  to  small  deviations  in  . In  fact, 
one  definition  of  the  entrance  length  is  that  point  at  which 
the  wall  friction  factor  reaches  a constant  value.  The 
maximum  values  of  turbulence  intensity  near  the  wall  mainly 
govern  wall  pressure  fluctuations  (p^^),  which  are  in  turn 
directly  related  to  the  wall  shear  stress  [19]: 

= PwV2.56  ' 

3 . Large  and  Fine  Scale  Turbulence 

Turbulence  is  a superposition  of  eddies  of  various 
sizes  and  energies;  there  are  no  distinct  separate  vortices. 

The  kinetic  energy  contained  by  these  various  sized  eddies 
is  determined  by  their  vdrticity  or  intensity  of  velocity 
fluctuations.  An  energy  spectrum  shows  how  the  energy  is 
distributed  according  to  the  various  frequencies.  In  general, 
eddy  size  decreases  as  the  wall  is  approached,  but  is  size- 
limited  by  the  fluid  viscosity.  Increased  velocity  decreases 
the  lower  size  limit,  while  the  upper  limit  is  determined 

I 

mainly  by  the  test  apparatus  size.  | 

An  inherent  problem  in  the  characterization  of  flow  j 

I 

situations  with  widely  varying  eddy  sizes  is  one  of  scaling. 

Ideally , by  including  certain  factors  in  the  empirical  results,  | 

one  should  be  able  to  apply  them  to  any  possible  flow  i 

situation.  But  considering  the  marked  differences  between  ! 

the  flow  within  a small  test  apparatus  and  along  the  side  of  ;! 
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a ship  (where  vortices  may  be  several  meters  wide),  it  is 
not  surprising  that  this  is  indeed  not  the  case.  However, 
controlled  small-scale  testing  is  still  of  great  value. 
Parameters  that  are  important  in  small-scale  tests  will  also 
be  important  in  service  situations,  although  not  to  the  same 
degree.  Knowledge  of  how  these  parameters  affect  corrosion 
is  of  great  value,  and  small-scale  testing  can  establish 
trends  and  reduce  the  requirement  for  full  scale  testing. 

The  energy  spectrum  shown  in  Figure  33  indicates  the 
range  of  eddy  size  (inverse  of  wavenumber  or  frequency)  and 
their  corresponding  energies.  Frequencies  in  a single  flow 
field  may  vary  from  1 to  10,000  cycles  per  second.  The 
eddies  that  contain  the  most  energy  are  neither  the  largest 
nor  the  smallest.  However,  the  larger,  more  permanent  ones 
may  still  contain  20  per  cent  of  the  total  energy.  As  eddy 
size  decreases,  dissipation  by  viscous  forces  increases,  up 
to  a maximum  for  a certain  sized  eddy  [l9].  Bremhorst  has 
shown  that  the  dominant  role  in  momentum  transfer  near  the 
wall  is  played  by  only  a limited  frequency  band  [58].  < In 
pipe  flow,  he  has  delineated  two  momentum  transfer  processes: 

(1)  A net  positive  transfer  in  the  high  frequency  range. 

(2)  A net  negative  tranfer  that  returns  low  momentum 
fluid  to  the  wall  at  the  lower  end  of  the  spectrum. 

Difficulties  arise  in  the  use  or  analysis  of  energy 
spectrums . Davies  has  pointed  out  that  a trace  is  not  like 
a fingerprint  for  turbulence;  it  doesn't  represent  even  a 
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Figure  33. 


Energy  spectrum  in  turbulent  flow. 
From  Hinze  [19]  . 
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blurry  one  [61]-  Even  in  one  specific  region,  a wide 
spectrum  will  be  present,  not  a single  peak  corresponding 
to  the  preferred  eddy  size  at  that  location.  Due  to 
averaging  that  is  done,  even  the  existence  of  a particular 
frequency  does  not  necessarily  imply  the  existence  of  a 
corresponding  identifiable  coherent  structure  in  the  flow. 

Finally,  even  an  array  of  organized  structures  of  similar 
scales  may  produce  a broad  spectrum. 

The  various  sized  eddies  may  move  at  different  velocities, 
and  be  inclined  to  the  wall,  some  from  5 to  15°  [19].  Mean 
shear  acts  to  amplify  preferred  eddy  structures  while  others 
decay . 

Large  eddies  interact  with  smaller  ones,  affecting  them 
in  several  ways: 

(1)  Deformation  of  the  smaller  by  shear  flow  of  the  larger; 

Moore  has  found  the  elliptical  deformation  of  vortices  due 
to  the  irrotational  straining  field  of  nearest  neighbors  [63]. 

(2)  Entrainment  of  the  small  eddies  by  the  flow  velocity 
caused  by  the  bigger  eddies,  particularly  for  eddies  of 
widely  differing  size. 

(3)  Transfer  of  energy  from  the  larger  to  the  smaller 
(accounted  for  by  viscosity)  , resulting  in  a cascade  process 
for  breakdown  of  eddies  . 

i 

The  structural  features  of  large  eddies  include: 

(1)  Pairs  of  contra-rotating  streamwise  vortices.  ' 
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(2)  Attached  eddies  elongated  and  flattened  in  the  stream- 
wise  direction,  with  centers  spaced  spanwise  approximately 
= 80  and  perpendicular  tothe  wall  y^  = 50  [19]. 

Naraynan  has  found  the.  width  of  the  fine  scale  region  to  be 
a function  of  Re;^  , either  sheetlike  vorticity  regions  with 
thickness  17  , or  filament-like  with  diameter  V and  spacing 
X [53]. 

These  large  eddies  can  be  considered  the  debris  of  ini- 
tially unstable  flow  [55].  They  are  known  to  play  an  impor- 
tant role  in  shear  layer  growth,  turbulent  mixing,  noise 
production,  and  turbulent  transport  of  momentum,  heat  or 
mass  [61]  . Wygnanski  and  Champagne  demonstrated  that  they 
possessed  an  ordered  structure  in  pipe  flow  and  grew  very 
slowly  if  at  all  [64].  Hinze  stated  that  the  diffusive 
action  of  turbulence  was  mainly  determined  by  the  larger 
eddies  [19].  Near  the  wall,  Lee  [65]  found  the  large  eddies 
to  be  greatly  elongated,  with  lateral  dimensions  on  the  order 
of  X^ / 2 , whereas  those  in  the  core  have  been  found  by  Hinze 
[19]  to  be  elongated  with  diameters  up  to  half  the  pipe 
radius  . 

Smaller  eddies  decay  faster  than  large  ones . They  have 
a larger  velocity  gradient  and  a larger  viscous  shear  stress 
counteracting  the  eddying  motion  [19].  In  flow  of  this  strong 
viscous  nature,  molecular  effects  dominate.  Ueda  demonstrated 
that  both  the  fine  and  large  scale  structures  in  the  laminar 
sublayer  are  influenced  by  the  inrush  phase  of  high  axial 
momentum  fluid  into  the  low  momemtum  fluid  [66]. 
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Lindgren  found  that  fine  turbulence  generated  at  the 
entry  decayed  with  flow  downstream,  until  a point  was  reached 
where  a coarser  equilibrium  structure  existed  [67].  Mulhearn 
observed  that  the  initial  structure  developed  by  the  generat- 
ing device,  such  as  a grid,  was  strongly  anisotropic.  However, 
downstream  the  structure  became  remarkably  similar,  differing 
only  in  the  length  scales,  dependent  on  the  particular  shear 
generating  device  [68]  . Thus  the  details  of  the  initial 
conditions,  i.e.  tripping  wire,  grid,  or  surface  roughness 
become  of  minor  importance,  except  in  scale  determination. 

Even  vortices  of  comparable  size  cannot  coexist  for 
long.  As  they  are  convected  downstream,  they  must  either 
disappear  or  amalgamate  [63].  Annihilation  of  neighboring 
vortices  occurs  as  they  tend  to  rotate  about  one  another. 
Roberts  and  Christian  found  this  to  occur  when  the  distance 
between  centers  reached  a value  less  than  1.7  times  the 
diameter  [69].  More  found  the  value  to  be  2.8,  at  which 
distance  mutual  distortion  caused  the  weaker  neighbor  to 
disintegrate  [63].  Its  vorticity  ingested  by  its  neighbors 
caused  vortices  to  grow  linearly  in  spacing  and  scale  with 
downstream  distance.  This  ingestion  of  irrotational  fluid 
would  require  turbulent  mixing  and  diffusion,  and  result 
in  a statistically  smooth  linearly  growing  mixing  layer  [55]. 

Hinze  has  found  that  the  irregularity  and  randomness 
of  real  turbulent  motion  are  the  main,  if  not  sole,  deter- 
minants of  the  transport  and  diffusion  of  turbulence  [19]. 
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Davies  has  emphasized  the  intermittency  produced  by  large 
structures  [6l],  while  Ueda  emphasized  that  the  inter- 
mittency of  the  fine  scale  structure  increased  with  Re;^ 

[66]. 

This  intermittency  is, however,  largely  absent  in  fully 
developed  pipe  flow,  making  the  burst  period  nearly  constant. 
The  repetitive  ordered  structures  do  strongly  interact  with 
each  other,  coalescing  as  they  move  with  the  flow.  These 
interactions  strongly  influence,  if  not  control,  many  of 
the  observable  effects  of  turbulence,  such  as  enhanced 
transport  process  . 

Hinze  [19]  has  shown  that  turbulence  in  the  high  wave- 
number  range,  i.  e.  smaller  eddies,  can  only  be  determined 
by  internal  parameters,  namely: 

(1)  Viscosity  v which  determines  the  rate  of  dissipation. 

(2)  Energy  flux  e , which  equals  the  rate  of  energy 
dissipation  per  unit  mass. 

Thus,  whereas  turbulence  in  the  core  is  affected  by  u 
and  the  curvature  of  the  velocity  profile  at  that  point, 
the  basic  parameters  in  the  wall  region  are  v and  u^. 

Studies  of  dissolution  patterns  on  a soluble  surface, 
like  plaster-of-Paris , give  additional  clues  to  the  effect 
of  hydrodynamics  on  corrosion  [70].  Just  as  with  sand  ripples 
on  the  ocean  floor,  as  water  flows  over  the  surface,  periodic 
troughs  or  flutes  develop.  Any  mass  transfer  boundary  layer 
that  may  exist  over  the  sample  surface  is  insufficient  in 
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itself  to  initiate  such  a periodic  pattern.  The  flow  itself, 
as  detailed  in  Figure  34,  acts  to  initiate  and  propogate  the 
contours.  Once  the  surface  irregularity  is  established,  the 
peaks  act  as  sites  for  flow  separation,  where  mass  transfer 
rates  are  low.  Flow  reattaches  downstream  where  mass 
transfer  rates  are  high.  Coupled  with  the  recirculation 
region  between,  troughs  are  formed  and  subsequently  spread 
downstream.  For  Reynolds  numbers  much  greater  than  that 
needed  for  flute  stability,  pot-holes  developed  in  the  slopes 
of  profiles,  where  flow  is  reattaching.  The  characteristic 
dimensions  of  the  resultant  dissolution  patterns  were  found 
to  depend  on  the  shear  stress,  fluid  viscosity  and  velocity, 
and  diffusivity.  The  analysis  is  complicated  by  the  fact 
that  the  dissolution  - rate  velocity  perpendicular  to  the 
wall  is  an  unknown  function  of  general  turbulent  flow. 

Experiments  by  Rudd  with  drag  reducing  polymers  further 
point  out  the  effect  of  fluid  variables  on  flow  properties 
like  shear  stress.  Adding  the  polymer  causes  earlier 
transition  to  turbulence,  thickens  and  stabilizes  the  vis- 
cous sublayer  (at  the  same  value  of  Reynolds  number),  but 
has  little  effect  on  turbulence  intensity.  Spanwise 
velocity  fluctuations  are  supressed,  accompanied  by  in- 
creased eddy  size  at  the  wall.  Only  10  to  100  ppm  can 
decrease  the  friction  factor  by  80  per  cent  [71]. 

Pipe  flow  would  certainly  seem  to  be  different  from 
that  along  a flat  plate.  Indeed,  the  pressure  gradient 
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Figure  34,  Hydrodynamic  regions  for  flow  over  one  flute. 

(1)  Main  turbulent  flow.  (2)  Laminar  free 
shear  layer.  (3)  Transition  to  turbulence 
and  entrainment.  (4)  Reattachment.  (5) 
Irregular  recirculating  flow,  (6)  Lee-slope 
boundary  layer.  (7)  Streamward-slope  bound- 
ary layer.  (8)  Lee-slope  boundary-layer 
separation  with  small  eddies  . From  Blumberp- 
[70]. 
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within  the  pipe  would  cause  the  flow  to  differ  significantly. 
Unlike  in  external  flow,  the  pressure  gradient  term  in  the 
Navier-Stokes  equation  is  non-zero.  This  prevents  direct 
comparison  of  the  two  flow  situations. 

In  a non-circular  duct,  secondary  motion  from  the  center 
towards  the  corners  of  only  1 per  cent  the  axial  velocity 
can  become  significant.  The  velocity  profile  is  bulged 
outward  near  the  corners,  as  Figure  35  shows.  Mailing 
showed  that  these  secondary  vortices  were  generated  by 
gradients  of  Reynolds  stresses  across  the  duct  cross  section 
[72].  Landolt  investigated  the  significance  of  these 
secondary  motions  on  high  rate  metal  dissolution  during 
electrochemical  machining.  His  results  using  a rectangular 
duct  were  consistent  with  those  in  a circular  duct,  despite 
the  fact  that  the  current  distribution  in  the  former  was 
inhomogeneous  [74].  This  lends  support  to  the  use  of  a 
non-circular  duct  for  corrosion  studies,  particularly  useful 
when  microscopic  examination  of  the  metal  surface  is  desired. 
However,  if  these  secondary  motions  have  a negligible  effect 
on  dissolution  rate,  this  would  seem  to  imply  that  small 
variations  in  turbulent  structure  would  be  similarly 
unimportant . 

Velocity  and/or  turbulence  affects  the  anodic  behavior 
of  a metal  in  two  ways : 

(1)  Removal  of  corrosion  products. 

(2)  Changes  in  concentration  of  chemical  species  at 
the  metal  surface . 
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The  manner  in  which  oxygen  is  broughc  to  the  surface,  i.e. 
to  the  local  cathodes,  must  be  considered.  It  seems  obvious 
intuitively  that  turbulence  parameters  either  affect  the 
corrosion  rate  in  themselves  or  affect  other  variables 
that  do.  That  the  corrosion  rate  varies  between  the 
laminar  and  turbulent  flow  regimes  in  circular  ducts  has 
been  pointed  out  by  Ross  [ 40]  . He  found  that  when  the  log 
of  corrosion  rate  was  plotted  against  the  log  of  Re^ 

(Reynolds  number  based  on  diameter)  two  linear  portions 
stood  out.  Where  Re^  was  less  than  2000,  the  corrosion  rate 
was  proportional  to  (Re^)^  , while  for  higher  values  it  was 
simply  proportional  to  Re^.  This  shows  that  a velocity 
change  in  the  laminar  regime  will  give  a smaller  change  in 
corrosion  rate  than  an  equal  velocity  change  in  the  turbulent 
regime.  However,  this  analysis  does  not  delineate  whether 
the  increased  corrosion  rate  is  due  to  the  increased  velocity 
(inherent  in  Re^  when  the  same  test  apparatus  is  used),  or 
to  some  tangible  or  intangible  in  the  turbulence  itself. 

To  accent  the  importance  of  turbulence  per  se  would 
require  experiments  on  separate  test  setups  where  the 
velocity  could  be  the  same,  but  Reynolds  numbers  vary. 

Varying  distances  from  the  leading  edge  of  a flat  plate 
would  seem  appropriate,  or  varying  the  cross  section  in  a 
flow  channel.  Cornet  similarly  found  the  increase  in  corro- 
sion rate  to  be  more  rapid  with  increasing  Reynolds  number  in 
the  turbulent  regime,  so  long  as  the  rate  was  controlled  by 
the  transport  of  oxygen  to  the  surface,  i.e.  diffusion 
controlled  [11] . 
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F.  HYDRODYNAMIC  CHARACTERIZATION 

In  the  present  work,  the  goal  was  to  compare  various 
turbulence  structures  or  patterns  and  connect  them  to  corro- 
sion phenomena.  Two  aspects  of  the  required  characteriza- 
tion are  the  size  differences  and  the  violence  of  the  motion. 
The  former  involved  a time  scale  and  a space  scale . The 
latter  is  the  turbulence  intensity,  u|  = (u^)^  for  the 
axial  component,  and  corresponding  quantities  for  the 
lateral  components . 

In  real  turbulence,  both  u^  and  the  velocity  gradient 
du^/dt  are  statistical  quantities,  so  rms  values  may  be 
used  as  averages.  As  the  wall  is  approached,  the  axial  ve- 
locity u^  decreases  while  its  fluctuation  u^  remains  nearly 
constant.  This  results  in  an  increase  in  the  turbulence 
intensity  u^/u^,  as  shown  in  Figure  36.  Yet  closer  to  the 
wall,  u^  continues  to  decrease,  with  u^  now  also  decreasing 
due  to  viscous  damping.  So  the  turbulence  intensity  decreases. 
A peak  then  results  at  some  characteristic  distance  from 
the  wall.  This  peak  region  is  usually  associated  with  the 
maximum  production  and  dissipation  of  turbulent  energy. 

Curves  found  by  various  investigators  are  shown  in  Figure  36 . 
The  peaks  are  seen  to  range  from  approximately  y*  ~ 12  to  17, 
with  maximum  values  of  u^  ranging  from  2.6  to  3.0.  This  is 
surprisingly  constant  over  quite  a range  of  Reynolds  numbers 
from  5600  to  35,500.  The  results  of  several  investigators 
on  the  limiting  value  of  turbulence  intensity  at  y+  = 0 range 
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from  0.24  (Eckelmann)  to  0.4  (Hinze),  with  Laufer  obtaining 
0.3  and  Mitchell  and  Hanratty  0.32.  While  the  use  of 
anemometer  probes  is  based  on  a basic  law  of  steady  state 
heat  transfer,  it  must  be  noted  that  real  turbulence  is 
anything  but  steady.  Perry  has  reported  that  differences 
of  25  per  cent  in  rms  turbulence  levels  is  not  uncommon  in 
the  literature  [75]. 

The  radial  component  of  the  fluctuating  velocity  gradient 
cannot  be  ignored,  certainly  when  considering  mass  transport 
to  the  wall.  Indeed,  the  fluctuations  in  the  direction  of 
mean  flow  have  little  effect  on  mass  transfer.  Sirkar  and 
Hanratty  showed  that  these  strong  transverse  velocities 
extend  right  up  to  the  wall  and  primarily  govern  the  fluc- 
tuating concentration  field  [76].  Hinze  has  found  that 
these  perturbations  in  velocity  that  propogate  towards  the 
wall  do  so  at  the  friction  velocity  u^.  With  the  sensing 
wire  perpendicular  to  the  flow,  both  the  axial  velocity 
component  u^  and  that  perpendicular  to  the  wall  U2  are 
being  measured. 

Although  the  turbulence  intensity  is  defined  as  u'/u, 
the  rms  voltmeter  reading  (e’)  from  the  hot  film  probe  does 
not  give  u'  directly.  It  can  be  calculated  if  a particular 
empirical  relationship  is  assumed  between  heat  transfer 
and  velocity.  Using  the  relationship: 
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gives : 
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2 . . 1/2 
where  e is  the  point  at  which  the  curve  of  u,  (deter- 
o 1 

2 

mined  from  the  rotometer)  versus  e (the  voltmeter  readout) 
crosses  the  e''  axis. 

A typical  curve  is  shown  in  Figure  37.  Readings  were 

taken  both  as  the  velocity  was  increased  and  again  while  it 

was  decreased.  Separate  runs  were  made  at  five  distances 

2 

from  the  wall.  The  resultant  values  of  e were  all  obtained 

o 

as  shown  in  Figure  37.  It  should  be  noted  that  the  open 
channel  gave  values  less  than  the  channel  with  a mesh  at 
the  entrance  for  the  two  positions  closest  to  the  wall  (0.23 
and  0.54  mm),  equal  values  at  1.09  mm  and  greater  for  the 
two  positions  farthest  from  the  wall  (3.0  and  1.8  mm).  In 
general,  the  turbulence  intensity  values  decreased  with 
increased  velocity.  At  a particular  velocity,  the  intensity 
increased  as  the  wall  was  approached,  as  shown  in  Table  6 
for  6 m/sec  (19.7  ft/sec).  The  values  did  not  peak  as 
expected  even  though  the  wall  was  approached  to  within  0.23 
' mm.  This  was  believed  due  to  the  disturbance  caused  by 

what  (at  that  small  a standoff  from  the  wall)  was  a relatively 
large  obstacle:  the  probe. 

Another  quantity  in  the  turbulent  flow  that  can  be 
derived  from  the  hot  film  anemometer  output  is  the  so-called 
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TABLE  6 

Distance  from  Wall  Turbulence  Intensity 


(mm) 

Open 

Channel 

Inlet 

Mesh 

3.0 

.030 

.033 

.030 

.038 

1.8 

.037 

.038 

.039 

.040 

1.1 

.047 

.047 

.047 

.048 

0.54 

.047 

.045 

.051 

.052 

0.23 

.050 

.048 

.054 

.051 

Values  of  turbulence  intensity  at  an  average  flow  velocity 
of  6 m/sec  in  an  open  flow  channel  and  with  an  inlet  mesh. 
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microscale  of  turbulence  ( X ) which  is  related  to  the  size 
of  the  smallest  eddies  in  the  turbulent  field,  those  eddies 
that  become  dissipated  due  to  viscosity  effects.  The 
signal  is  processed  through  a differentiator  with  a time 
constant  t and  read  out  on  an  rms  voltmeter.  This  value  (r) 
in  units  of  velocity  is  then  used  in  the  expression: 


u 


r 


As  with  the  turbulence  intensity,  the  mesh  did  have 
an  effect  on  r.  Except  for  the  position  closest  to  the 
wall,  the  mesh  readings  were  consistently  higher  than 
those  taken  in  the  open  channel.  This  would  decrease  the 
value  of  X , exactly  the  effect  expected.  The  character 
of  any  fluid  dynamic  process  like  turbulence  is  directly 
determined  by  the  ratio  between  inertial  and  viscous  forces . 
Therefore,  a Reynolds  number  of  turbulence  can  be  formed 
using  this  microscale  as  the  characteristic  length  and 
u or  u ' as  the  velocity  term. 

G.  MASS  TRANSFER 

Heat  and  mass  transfer  are  analagous  processes.  In 

any  boundary  layer  dominated  process,  it  has  been  experi- 
l/3 

mentally  determined  that  Nu  a Sc  . The  Sherwood  number 

(Sh)  is  the  mass  transfer  equivalent  of  the  Nusselt  number. 
1/3 

Plotting  Sh/Sc  versus  Re^  gave  Tvarusko  approximately 
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linear  results  with  horizontal  cylinders  over  a wide  range 
of  Re^  up  to  4 X 10^  [77].  He  found  an  increase  in  the  ratio 
of  only  15  per  cent  as  a result  of  increasing  the  turbulence 
intensity  from  0.15  to  1.0,  and  this  much  change  occurring 
only  at  the  highest  Reynolds  numbers.  Sirkar  has  shown  that 
at  high  Schmidt  numbers,  the  mass  transfer  is  actually 
controlled  by  this  fluctuating  velocity  field  [76].  The 
exponent  of  the  Schmidt  number  may  vary,  depending  on  the 
relative  importance  given  to  the  transfer  of  mass  by 
momentum  or  diffusion.  Leumer  has  summarized  the  dimension- 
less mass  transfer  correlations  most  prevelant  in  the 
literature  [43].  The  thickness  of  the  mass  transfer  boundary 
layer  ( 8^)  stabilizes  very  rapidly  with  distance  along  a 
wall.  In  the  present  flow  channel,  for  velocities  of  4 to 
8 m/sec,  it  is  on  the  order  of  20Mm  thick,  based  on  the 
relationship  given  by  Wranglen  [78]. 

The  variation  in  boundary  layer  thickness  does  correlate 
well  with  certain  parameters  that  have  a known  effect  upon 
corrosion.  A thinner  boundary  layer  has  a steeper  concentra- 
tion gradient.  However,  to  assume  that  this  always  results 
in  increased  corrosion  would  be  mistaken.  Its  thickness 
can  be  altered  in  two  different  ways: 

(1)  Increasing  the  velocity,  i.e.  Reynolds  number, 

decreases  8_  and  does  increase  the  corrosion  rate, 
m 

(2)  Increasing  the  diffusion  coefficient  (by  an 

2 

increase  in  temperature  for  example,  which  also  decreases  t'  ) 
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decreases  the  Schmidt  number  ( t>  /D,-,  ) and  therefore  increases 

-1/3 

the  boundary  layer  thickness.  But  since  a Sc  , Dg 

increases  at  a faster  rate  than  8 . The  mass  flux  is  pro- 
portional  to  Dg  / 8^^,  so  the  net  effect  is  an  increase  in 
corrosion  rate. 

Recapitulating,  in  (1)  above,  a decrease  in  5^  leads  to 

an  increase  in  corrosion  rate,  while  in  (2)  an  increase  in 

8 leads  to  an  increase  in  corrosion  rate, 
m 
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V . SUMMARY 


As  a result  of  this  research,  the  following  conclusions 
concerning  90/10  Cu-Ni  have  been  reached: 

(1)  A flow  channel,  as  opposed  to  other  test  setups, 
has  the  important  advantages  of  allowing  the  specimens  to 
be  easily  instrumented  for  electrochemical  monitoring, 
ensuring  fully  developed  turbulent  flow  across  the  specimen, 
and  permitting  continuous  viewing  of  the  specimen  surface. 

(2)  When  crevices  or  steps  perpendicular  to  the  flow 
are  avoided,  test  runs  for  48  hours  show  no  preferential 
dissolution  at  the  leading  edge  that  might  be  due  to  the 
establishment  of  the  mass  transfer  boundary  layer. 

(3)  There  is  a strong  dependence  of  the  corrosion 
rate  of  90/10  Cu-Ni  on  temperature,  attributed  mainly  to 
increased  reaction  kinetics  and  the  increase  in  the  diffusion 
coefficient  of  oxygen. 

(4)  The  corrosion  rate  of  90/10  Cu-Ni  increases  non- 
linearly  with  ve^-ocity,  ranging  from  6.5  mpy  at  4 m/sec  to 
13  mpy  at  7.8  m/sec.  (All  the  corrosion  rates  are  expected 
to  be  on  the  high  side,  due  to  averaging  over  the  relatively 
short  duration  of  48  hours,  when  the  corrosion  rate  is 
expected  to  be  higher  than  steady  state.) 

(5)  The  turbulence  intensity  was  investigated  using 
hot  film  anemometry.  It  increased  as  the  wall  was  approached 
down  to  a minimum  distance  of  0.23  mm.  At  that  minimum 
distance  tested,  the  values  ranged  from  .046  to  .085,  roughly 
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doubling  as  the  flow  velocity  was  decreased  from  7 m/sec 
to  2.5  m/sec.  Introducing  a mesh  at  the  channel  entrance 
(51  hydraulic  diameters  away  from  the  specimen)  increased 
the  turbulence  intensity  roughly  2 per  cent,  but  at  the 
only  velocity  tested  (4  m/sec)  had  no  effect  on  corrosion 
rate  . 

(6)  Empirical  scaling  factors  are  necessary  to  relate 
the  variable  effects  of  hydrodynamic  parameters  on  corrosion 
for  widely  varying  flow  situations.  These  are  difficult  to 
develop  experimentally,  but  careful  modeling  of  the  rela- 
tionship in  a small  test  facility  can  at  least  establish  the 
trends  in  a service  situation. 

(7)  Fluctuating  local  peaks  in  wall  shear  stress  may  be 
orders  of  magnitude  greater  than  the  average  shear  stress 
calculated  from  macroscopic  hydrodynamic  conditions . This 
fluctuating  stress  may  affect  corrosion  in  several  ways, 
such  as  altering  the  form  and  distribution  of  anodic  films, 
limiting  the  buildup  of  corrosion  products  and  removing 
existing  accumulations  of  corrosion  products. 

(8)  The  concentration  of  electrochemical  products  (metal 
ions)  and  reactants  (oxygen)  at  the  metal/electrolyte  inter- 
face is  controlled  by  molecular  diffusive,  natural  convective 
and  eddy  diffusive  (forced  convective)  modes  of  mass  transport. 
The  model  for  corrosion  in  a flow  field  must  particularly 
include  the  effect  of  eddy  diffusivity,  which  can  often  be 

the  major  component  of  mass  transport.  When  considering 
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eddy  diffusivity,  the  details  of  the  flow  structure,  in 
terms  of  the  eddy  size  distribution  in  the  turbulent  flow 
field,  must  be  taken  into  account,  as  this  is  a real 
physical  difference  in  various  flow  situations . The  model 
by  which  this  feature  is  accounted  for  might  incorporate 
the  effects  of  mixing  in  dimensionless  groups  like  the 
Sherwood  and  Schmidt  numbers,  or  might  adapt  experimentally 
measurable  parameters  that  characterize  the  turbulence, 
such  as  turbulence  intensity  and  microscale. 
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